y,fllecase 

P 5febl972, 

file blocks 

punch 
parity 
routines 3 

text print 






aM 



\^ l < _r. V 



3' 



^ 



iara-<3 am 


1 


nam 


1 


decimal 




Input 




fancy 


1 


simple 


1 


newfloat 




pack 


12 


sqrt 


1 


log 


2 


exp 


2 


cossln 


2 


elliptic 


3 


arc sine os 


2 


at an 


2 


cbrt 


1 


ran 


1 


phi 


1 


afb 


1 


bessel 


6 


gauss 


h 


afl 


l 


gamma 


2 


descrlptl 


io 

— s j- 


— atuie plUx; — 


J 


phlnvr 


2 


shellsort 


2 


dlsplSmar 


3 


ktmdesc 


4 




T. 


ktra 


6 


plotter 




minimal 


1 


text .plot 


11 


line 




shorter 


2 


basic 


3 


dpy 


3 


ascil 




flexo 




program 




b 


3 


a 


11 


subroutln 


h 


prime 




tester 


l 


pen 




follow 




old 


l 


new 


2 


random 




number 




tester 





ms 


2 


earrelati 


2 


simple 


1 


more 




random 


1 


gausslan 


3 


single 




precision 




sQrt 


1 


sQLrt.sum 




sQrs 


1 


sinVcos 


1 


sine os 


1 


lognat 


1 


double 




precision 




sqrt 


1 


decimal 




print 




fraction 




old 


4 


integer 


1 


print 




octal 


1 


decimal 




signed 


1 


scope 


»—-. 


decimal 




print 


2_ 


line 




fast 


1 


window 




fixed 


2 


variable 


3 


character 




display 


5 
6 


faster 


7-bit 


2 


et 




typewrite 




simulator 


3 


read /writ 


7 



ptAi/\"{~ 






505 free tape blocks 

free directory entries 



Numeric type out routines 
5 February 197 2 



1.1 

Numeric type out routines 



/octal print subroutine 

/leading zero suppressing 

/call by jdp opt with number in AC 

opt, o -U>cW- ft- 

lio opa+4 ax*-****M(. ^ <ji 

c -^rc"l 3s^ > \ . | 

sni ^ LKC S 17l>&^ fwecU^ 

Jmp . --2 
opa, c^qATp >J— xA~ 

law mil i i ii m " 
ivk 100 
opa+4, law 

C ^rcT3s*-> — _ . p 

jmp i opt 



Numeric type out routines 

/ decimal print subroutine 

/ called by jdp dpt with number in AC 






eNi£^> « r X ^S. J^^acfC 



dpt, 



dpp, 




lac^^T 
dzm dp2 
dac dp3 
mul (1 

div .+1 
12 

sas dp2 
jmp dpp 
T3AP| 
law 
ivk 100 
lac dp3 
dac dp2 

lac dpi 
sas dp2 

jmp dpp 
jmp 1 dpt 



H 



tf AA?#P 
{Ja A( 



variables 
constant 



Decimal input routines 
5 February 197 2 



J. « JL 

Decimal input routines 

/unsigned decimal input from typewriter 

/call with <jdp gnu> 

/returns with number in A, non-numeric break character in I, X, 

/rph 4aug7^ 

gnu, 

Z IX P f 

gnvj mul (5 -•*■** ,0 * t.fe* ^^ '' ' K 

law i 20 

XMX "\ 

X+IX x Uui 2.1 

TAI 

ivk 200 

X;+I<= <fc vnake VUiS X +-X*^ 

TXXlf 
jmp gnv 
jmp i gnu 



2.1 
Dec lira 1 input routines 

/signed decimal input from typewriter 

/call with <jdp decin> 

/returns with number in A, non-numeric break character in i, x. 

/restarts on centerdot, accepts minus sign before number 

/rph 4aug7^ 

decin, — 

decin+lj dzm gny /!- 

ZIAXP 'J?Y' % c r *' |4 ~ 

gntj lac gnz 

A~II_ 

dio gnz 
gnvj mul T5 

law i 20 

XAAX 

X+EC 

TAI 

ivk 200 
sas (40 
jmp gnw 
law 77 
ivk 100 
jmp decin+1 
gnw 5 sas \y\ 

jmp gnx 

lio gny 
CUM] 

jmp gnt 
gnx j dio gny 

XhAXP 
X+I<= 
IXXlj 

jmp gnv 

xor gnz 
jmp i dec in 



/knob hysteresis 

/iam or dam 

/elm or e2m 

/call with 

/ ckn nOO ^ 

/ r Jsp hys -^ 

-/kbn, (^ ' 

/expects old value of knob in kbn, 
/ returns new value there, and in AC 

hys, TAX 

lac i / "5<d- J<k ucJu*. 

AMIAI< / s^ **v 

NAA 

add (2 /increase for more hysteresis 

SAA< 

ZAl 

adm i 0^ \ ^ /i^ !**WWX e*A U JteUled \C «* **W< u 



Text printing routines 
5 February I972 



1.1 

Text printing routines 

/text printing subroutine 

/call by jsp txx followed by text 

/text should end with • 

/control returns to location after end of text 

/iam or dam 

/runs in core only 

txx, dap txy 

aarn 

11 o txy 

idx txy 

lac (607600 

rcl 6s 

sad (lse£- ~?c, 00 d , 
txy, jmp • 

sad (swp 

jmp txx+1 

ivk 100 

jmp txy-3 



2.1 
Text printing routines 

/text printing subroutine 

/call by jsp txx followed by text 

/text should end with • 

/control returns to location after end of text 

/nam 

/runs in core only 
txx, 



txy, 



dap 


txy 






lio 


i txy 






idx 


txy 






lac 


(6076OO 




rcl 


6s 






sad 


(lai 


-76 


viS-tO 


jmp 


• 






sad 


(swp 


*?fr 


troUO 


jmp 


txx+1 






ivk 


100 






jmp txy-3 







7-bit character routines 
5 February I972 



1.1 



7-bit character routines 



/7-bit character put and get 

/5 seven-bit characters packed in a word pair 

/bit-0,word spare 

/chr in word (hits 1-7) 

/chr l in word ^bits 8-l4 ) 

/chr 2 in word ibits 15-^-7) a ^d word a. ^bits 0-3 ) 

/chr J in word x ^ bits 4-ilo; 

/chr 4 in word 1 (bits II-I7) 



/get 

/character address in A, less than 400000 octal 

/character returned in A 
/&3 cycles 



get, 



/put 




elf 1 

jmp pgl 



/character address in A, less than 400000 octal 
/character in to 
/uses tl,t2 
/32 cycles 



put, 
Pg 1 , 



Pg 2 > 



Pg3, 



sht, 



ust, 








stf 


1 


mul 


(146315 


dac 


tl 


AMU 


TAX 




rcl 3s 


sza 


i 


law 


7 


lio 


i 1 


Ixr 


i 


X-*A2 




xct 


i sht 


szf 


i 1 


jmp 


Pg3 


and 


(-I77 


ior 


to 


CXX 




xct 


i ust 


Ixr 


tl 


dac 


i 


dio 


i 1 


jmp 


1 put 

(177 


and 


jmp 


i get 


rar 


8s 


ral 


7s 


rcl 4s 


ral 3s 


nop 




rar 


3s 


rcr 


4s 


rar 


7s 



/word address 
/correction factor 
/I [0-2] = 7,4,1,6,3 



/word 1 of pair 
/word of pair 

/shift into A [11-17] 

/get 

/leave others 

/put in the character 

/read ust backwards 

/unshift it 

/put pair back 
/mask it out 



/chr 0, unshift 

/chr 3, unshift 

/chr 2, shift 

/chr 1, unshift 

/chr 4, shift and unshift 

/chr 1, shift 

/chr 2 , unsh if t 

/chr 3 , shift 



1.2 
7-bit character routines 

ral 8s /chr 0, shift 



Ascii to f lexo conversion 
5 February 1972 



1.1 

Ascli to f lexo conversion 

ascii-et 

/ type n to not convert lower and upper case 

/ type c to undo n 

/ type r f or ascii reader to f lexo text 

/ type p f or f lexo text to a sell punch 



0/ 
100/ 



jmp ice 



jmp 102 
nop 
law 77 

ivk 
jmp In 



ftoa, 

sas (72 



sad (74 

jmp f cs 
ior cas 
TAX 

law 177 

and i tab 
sad (15 
jmp cr 
eve, skp 600 /skip if converting to upper case 

jmp i ftoa 

sub (14 

spa 

add (40 

add (100 

jmp i ftoa 
cr , aam 

xct ftoa 

law 12 /generate line feed after carriage return 

jmp i ftoa 
fes, sub (72 /case shift 

ral 5s 

dac cas 

idx ftoa 

jmp i ftoa 



2.1 



Ascii to f lexo conversion 



at of , 



af3, 



afl, 



af2, 
IT, 





sad (12 
jmp If 
elf 6 
sad (15 
stf 6 
TAAX 
xct eve 
jmp afl 
sub (100 
TAA> 
jmp afl 
sub (32 
TA> 
law 40 
A+XX 

lac i tab 
ell 
sar gs 
dac ftemp 
and (100 
sad cas 
jmp af2 
dac cas 
rar 5s 
add (72 
aam 

xct at of 
law 77 
and f temp 
jmp 1 at of 
szf i 6 
jmp af 3 
elf 6 
idx at of 
jmp 1 at of 



/last character was carriage return 
/skip if converting to lower case 



3.1 
Ascii to f lexo conversion 



define z f ,a 




fxlOOO a 


te rmin 




tab, z 


76,4o 


z 


76,61 


z 


76,62 


z 


14,63 


z 


76,64 


z 


76,65 


z 


76,66 


z 


76,67 


z 


75,70 


z 


36,71 


z 


77,0, 


z 


76,14 


z 


13,3 


z 


77,0 


z 


34,0 


z 


35,0 


z 


76,60 


z 


76,57 


z 


76,163 


z 


76,164 


z 


76,165 


z 


76,166 


z 


76,167 


z 


76,170 


z 


76,171 


z 


76,172 


z 


76,0 


z 


76,54 


z 


76,17 


z 


76,16 


z 


76,11 


z 


76,0 


z 


0,72 


z 


105,152 


z 


101,153 


z 


103,154 


z 


103,155 


z 


104,156 


z 


106,157 


z 


102,160 


z 


57,161 


z 


55,162 


z 


173,0 


z 


154,0 


z 


33,55 


z 


54,51 


z 


73,73 


z 


21,50 


z 


20,0 


z 


1,141 


z 


2,142 


z 

z 


3,143 
4,144 


z 


5,145 


z 


6,146 


z 


7,147 



3.2 

Ascii to f lexo conversion 

z 10,150 

z 11,151 

z 40,0 

z 56,56 

z 10*7,0 

z 133,10 

z 110,0 

z 120,40 /100 

z 161,42 

z 162,47 

z 163,44 

z 164,45 

z 165,41 

z 166,46 

z 167,74 

z 170,76 

z 171,136 

z 141,0 

z 142,0 

z 143,0 

z 144,0 

z 145,0 

z 146,0 

z 147,100 

z 150,77 

z 151,123 

z 122,124 

z 123,125 

z 124,126 

z 125,127 

z 126,130 

z 127,131 

z 130,132 

z 131,0 

z 157,75 

z 156,17 

z 155,16 

z 111,11 

z 140,0 

z 102,137 

z 61,112 

z 62,113 

z 63,H4 

z 64,115 

z 65,H6 

z 66,117 

z 67,120 

z 70,121 

z 71,122 

z 41,0 

z 42,0 

z 43,53 

z 44,135 

z 45,134 

z 46,133 

z 47,0 

z 50,101 

z 51,102 

z 22,103 



3.3 
Ascii to f lexo conversion 

z 23,104 
z 24,105 
z 25,106 
z 26,107 
z 27,110 
z 30,111 
z 31,0 
z 157,52 
z 156,0 
z 155AO 
z 103,0 
z 76,15 



4.1 
Ascii to f lexo conversion 

/call with jdp rset, AC contains field number of text 

/ac,io,f lags, address mode saved. Enter in iam,nam,dam 

/rbl must be a submultiple of 400 and a multiple of 40 

rbl=400 

dimension rbf (rbl) 

rset, 

dap rsr 

dac ac 

dio io 

rar 6s 

add (400-rbl 

dac ddp 

dzm rpc 

lac (isp rpc 

dac rfg /flag for end of file 

law rbf-frbl-1 /last word of buffer 

dap inp 

law rbf 

mta 

lio (340 

law 40 
rsr, ivk • 

hit 

lac ac 

rar 6s 

add rbf '& 

sub (20000 

dac dne 

lio io 

lac ac 

jmp i rset 



5.1 
Ascil to f lexo conversion 



jmp rp2 after end of file 



readch, 







rfg* 


isp 


rpc /becomes 




jmp 


inp 




law 


i 3 




dac 


rpc 




idx 


inp 




dap 


inq 




sas 


(lio rbf+rbl 




jmp 


9k 




law 


rbl 




a dm 


ddp 




lia 






ral 


6s 




dap 


rdr 




law 


rbf 




dap 


inp 




dap 


inq 




mta 






law 


rbl 


rdr, 


ivk 
hit 
ZAP 




9K 


sub 


(lio rbf 




add 


ddp 




sad 


dne 




jmp 


rp3 


inp, 


lio 
cla 






rcl 


6s 


inq, 


dio 






sas 


(77 




sad 


(13 




dzm 


rpc 


rpl, 


lia 






jmp 


i readch 


rp3, 


lac 


(jmp rp2 




dac 


rfg 


rp2, 


law 


14 




jmp 


rpl 



6.1 
Ascii to f lexo conversion 

/et buffer write 

/of 1 must be a submultiple of 4 00 and a multiple of 40 

/call with jdp wset, ac contains field number of text to be written 

of 1=400 

dimension ofb(ofl) 

wset, 



dac 


ac 


dap 


wrdr 


mta 


300 


nop 




sni 




bpt 




rar 


6s 


dac 


wf Id ~ 


add 


_(400-of 1 


dac 


odp 


law 


i 3 


dac 


cp 


law 


ofb 


dap 


ofp 


dap 


ofq 


lac 


ac 


dzm 


cist"" 


jmp 


i wset 



/last character 



7.1 
Ascii to f lexo conversion 

/call with jdp writec, io contains character in ivk format 

wlOO, /write from ac 

dac ac 
dio 95P 
jdp writec 
lac ac 
lio 95p 

jmp i wlOO 

w300, 

dac 95P 

dio ac 

lai 

jdp writec 

lac 95P 

lio ac 

jmp i w300 

writec, 

and i.77 

dac 9op 

sas (13 

jmp . 4 

law 77 /make sure this is preceded gy 77 

sas cist 

jmp force 

lac (lio ofb+of 1 

sas ofp 

jmp ofp-1 

law of 1 /write out buffer 

adm odp 

lia 

and (770000 

sas odp 

jmp wo5 

ral 6s 

mta 300 

nop 

sni 

bpt 

jmp . +2 
wo5> ral 6s 

dap wrd 

law of Id 

dap ofp 

dap ofq 

mta 

law of 1 20 
wrd, ivk 

hit 
ofp-1, lac 96p 

ofp, lio 

rar 6s 

rcl 6s 
ofq, dio 

isp cp 
jmp ou3 
law i 3 



7.2 

Ascii to f lexo conversion 

dac cp 

idx ofp 

dap ofq 

wox, lac 96" P 

dac cist /character just written 

jmp 1 writec 

ou3, lac 96p 

sas (77 



sad (13 
jmp ofp 
jmp wox 

force, lac writec 

dac acw 
law 77 
jdp writec 
lac acw 
dac writec 
lac ac 
jmp writec 1 



8.1 
Ascii to f lexo conversion 

/call with jdp clean, finishes up output, writes out bufer 

clean, _ 

dac acl 

dio io 

law 13 

sas cist /force 13 to precede end of text 

jdp wlOO 

lac ofp 

add odp 

sub wf Id 

add (20000+of l-[lio ofb] 

dac b?7 

law ofb+ofl 

dap ofp 

jdp writec /write out last buffer 

law b77-32 

mta 

lio (340 

law 60 
wrdr, ivk • 

hit 

lac acl 

lio io 

jmp i clean 

b77, 

20400 



9.1 
Ascii to f lexo conversion 

in, iam 

elf 7 

ivk 2 00 

sad (char c 

jmp con 

sad (char n 

jmp con-1 

sad (char rr 

jmp red 

sad (char p 

jmp pun 

jsp txx 

text ..35.?. 347740.. 

jmp in 

red, lac (30017 

mta 306 

jmp bsy 

law 2 

jdp wset 
11, ivk 17 

jmp re dun 

stf 1 

sza 

sad (377 

jmp il 

and (177 

jdp at of 

jdp wlOO 

jmp il 
re dun, szf i 1 

jmp roff 

jdp clean 
end, law 17 

mta 204 

dsm 
bsy, jsp txx 

text ..3577. busy. 34 7740.. 

jmp in 
roff, jsp txx 

text ..77-turn on reader and try again. 7740.. 

law 17 

mta 204 

jmp in 

con-1, ZAP 

con, law 600 

dap eve 

jmp in 

pun, 3a c (40017 

mta 306 

jmp bsy 

law 2 

jdp rset 

law i 20 

jdp fee 
ol, jdp readch 

sad (14 

jmp oe 



9.2 

Ascii to f lexo conversion 





jdp f toa 




jdp out 




jmp ol 


oe, 


law i 60 




jdp fee 




jmp end 


fee, 







TAI 




ivk 17 




SUP 




jmp .-2 




jmp i fee 


out, 







TAI 




law 2325 




rcr 7s 




ral 7s 




dap .+1 




ril 




rar 7s 




rcl 7s 




law 200 




A~IA 




ivk 17 




jmp i out 


txx, 


dap txy 




aam 




lio txy 




idx txy 




lac (6076OO 
rcl 6s 






sad (lai 


txy, 


jmp . 




sad (swp 




jmp txx+1 




ivk 100 




jmp txy-3 


constants 




variables 




start in 





10.1 
Ascii to f lexo conversion 



/ascii to f lexo / f lexo to ascii 

/runs in iam or dam 

/calling sequence for f lexo to ascii 

/(ascii to flexo is similar) 

/f lexo character in A 
/ jdp f tea 

/ jdp put 

/return here 

/"put 11 is a user supplied routine which accepts the 
/ ascii characters produced by ftoa in A, 

/ It may he called 0, 1, or 2 times 

/ during a single call to ftoa 

/these routines do not save any registers 

/"cas tt will contain the flexo case (0 or 100 ) 

/*cvc u should he a n skp n if conversion between upper 

/ case ascii and lower case flexo is not desired. 

/ It should be a l *skp 600 1 * if this conversion 

/ is desired, (this is the default condition) 

/at of uses flag 6 (it should be initially cleared) 



11.1 

Ascii to f lexo conversion 



ftoa, 



cvc, skp 600 /skip if converting to upper case 



cr 








sas 


g 


sad 


jmp 


fcs 


ior 


cas 


TAX 




lav; 


177 


and 


i tab 


sad 


(15 


jmp 


cr 


skp 


600 


jmp 


i ftoa 


sub 


(140 


spa 




add 


(4o 
(loo 


add 


jmp 


i ftoa 


aam 




xct 


ftoa 


law 


12 


jmp 


i ftoa 
(72 


sub 


ral 


5s 


dac 


cas 


idx 


ftoa 


jmp 


i ftoa 



/generate line feed after carriage return 
fcs, sub (7 s /case shift 



Ascii to flexo conversion 






atof, 



af3, 



afl, 



af2, 
If, 





sad (12 

jmp If 

elf 6 

sad (15 
stf 6 

TAAX 
xct eve 
jmp afl 
sub (100 
TAA> 
jmp afl 
sub (32 
TA> 
law 40 
A+XX 

lac i tab 
ell 

sar gs 
dac ft e rap 
and (100 
sad cas 
jmp af 2 
dac cas 
rar 5s 
add (72 
aam 

xct atof 
law 77 
and f temp 
jmp 1 atof 
szf 1 6 
jmp af 3 
elf 6 
idx atof 
jmp 1 atof 



/last character was carriage return 
/skip if converting to lower case 



13.1 
Ascii to f lexo conversion 



define z 


f ,a 






fxlOOO a 


termin 






tab, 


z 


76,4o 




z 


76,61 




z 


76,62 




z 


14,63 




z 


76,64 




z 


76,b5 




z 


76,66 




z 


7b, 67 




z 


75,70 




z 


36,71 




z 


77, o, 




z 


76,14 




z 


13,3 




z 


77,0 




z 


34,0 




z 


35,0 




z 


76,60 




z 


7§*57 




z 


76,163 




z 


76,164 




z 


76,165 




z 


76,166 




z 


76,167 




z 


76,170 




z 


76,171 




z 


76,172 




z 


76,0, 




z 


76,54 




z 


76,17 




z 


76,16 




z 


76,11 




z 


76,0 




z 


0,72 




z 


105,152 




z 


101,153 




z 


103 ,154 




z 


103,155 




z 


104,156 




z 


106,157 




z 


102,160 




z 


57,161 




z 


55,162 




z 


173,0 




z 


154,0 




z 


33,55 




z 


54,51 




z 


73,73 




z 


21,50 




z 


20,0 




z 


1,141 




z 


2,142 




z 


3,143 




z 


*T • *U*"T*T 




z 


5,145 




z 


6,146 




z 


7,147 



Ascii to f lexo conversion 

z 10,150 
z 11,151 
z 40,0 
z 56,56 
z 107,0 
z 133,10 
z 110,0 
z 121,15 
z 120,40 
z 161,42 
z 162,47 
z 163,44 
z 164,45 
z 165,41 
z 166,46 

z 167,74 
z 170,76 
z 171,136 
z 141,0 
z 142,0 
z 143,0 
z 144,0 
z 145,0 
z 146,0 
z 147,100 
z 150,77 

z 151,123 
z 122,124 
z 123,125 
z 124,126 
z 125,127 
z 126,130 
z 127,131 
z 130,132 
z 131,0 

z 157,75 
z 156,17 
z 155,16 

z 111,11 

z 140,0 
z 102,137 
z 61,112 
z 62,113 
z 63,114 
z 64,115 
z 65,116 
z 66,117 
z 67,120 
z 70,121 

z 71,122 
z 41,0 
z 42,o 
z 43,53 
z 44,135 
z 45,134 
z 46,133 
z 47,0 
z 50,101 
z 51,102 
z 22,103 



13.2 



/100 



13.3 
Ascii to f lexo conversion 

z 23,104 

z 24,105 

z 25,106 

z 26,107 

z 27,110 

z 30,111 

z 31,0 

z 157,52 

z 156,0 

z 155AO 

z 103,0 

z 76,15 



ET reading and writing routines 
5 February I972 






ET reading and writing routines 
typewriter to e.t. text 



dimension ac (3 ) 
setup, 



io=ac+l 



finish, 
eno, 



enn 



ent, 



epc 



var 



vars=. -1 



lio (ent 
law 
mta 307 
bpt 

dap enn 
ral 6s 
lia 

mta 405 
bpt 

dap eno 
mta 4 01 
law 7 
jdp wset 
jmp i setup 



law 

mta 4 01 
mta 204 
law 

mta 204 
jdp clean 
jmp i finish 

dap epc 
law 1700 
AAIA 
rar 6s 
dap var 
lio (ac 
law 11 
ivk 
law 1600 

and ac 
lia 

law 
TAXP| 
SIX 
iam 

xct i vars 
jdp writec 
law 71 

xct epc 
qit 



lac ac 
hit 
lac io 
repeat 14, hit 



/100 - out from ac 
/300 - out from io 



2.1 
EI} reading and writing routines 

/et buffer write 

/of 1 must be a submultiple of 4 00 and a multiple of 40 

/call with jdp wset, ac contains field number of text to be written 

of 1=400 

dimension ofb(ofl) 

wset, 

dac ac 

dap wrdr 

mta 300 

bpt 

rar 6s __ 

dac wf Id 

add £400-of 1 

dac odp 

law i_ 3 

dac cp 

law ofb 

dap ofp 

dap ofq 

lac ac __ 

dzm cist /last character 

jmp i wset 



3.1 
ET reading and writing routines 



writec, 









and 


i77 




dac 


96p 




sas 


(13 




jmp 


. 4 




law 


77 




sas 


cist 




jmp 


force 




lac 


(lio ofb+ofl 




sas 


ofp 




jmp 


ofp-1 




law 


ofl 




a dm 


odp 




lia 






and 


(770000 




sas 


odp 




jmp 


wo5 




ral 


6s 




mta 


300 




bpt 






jmp 


. 42 


wo5, 


ral 


6s 




dap 


wrd 




law 


ofb 




dap 


ofp 




dap 


ofq 


\ 


mta 






law 


ofl 20 


Wrd, 


ivk 
hit 




ofp-1, 


lac 


96p 


ofp, 


lio 






rar 


6s 




rcl 


6s 


Ofq, 


dio 






isp 


cp 




jmp 


ou3 




law 


i 3 




dac 


cp 




idx 


ofp 




dap 


ofq 


WCK, 


lac 


96p 




dac 


cist 




jmp 


i writec 


ou3, 


lac 


?6p 




sas 


77 
(13 




sad 




jmp 


ofp 




jmp 


WOK 


force, 


lac 


writec 




dac 


acw""" 




law 


77 




jdp 


writec 




lac 


acw 




dac 


writec 




lac 


ac 




jmp 


writec 1 



/make sure this is preceded "by 77 



/write out buffer 



/character just written 



4.1 

ET reading and writing routines 

/call with jdp clean, finishes up output, writes out bufer 

clean, _ 

dac acl*~ 
dio io 
law 13 

sas cist /force 13 to precede end of text 

jdp writec 
lac ofp 
add odp 
sub wf Id 
add (20000+of l-[lio ofb] 

dac b77 

law ofb+ofl 
dap ofp 

jdp writec /write out last buffer 

law b?7-32 
mta 

lio (340 
law 60 
wrdr, ivk . 

hit 

lac acl 
lio io 
jmp i clean 

20400 



5.1 

ET reading and writing routines 

constants 
variables 
start 



6.1 
ET reading and writing routines 

/e.t. text reader /writer 

/this is a set of subroutines to allow one to read and write e.t. type 
/text fields on drum using subroutine calls that work in the same way 
/as typewriter ivk instructions, i.e., they take and return characters 
/in the rightmost byte of the ac or io registers. These routines 
/preserve the contents of the io,ac,xr, and flag registers. 

/for writing onto the drumi 

/ jdp wset — this initializes everything for writing on the drum 

/ the ac must contain the field number on which text 

/ is to start. This may be any field. 

/ jdp wlOO — this takes one character from the ac and 

/ writes it on the drum 

/ jdp w300 — this takes one character from the io and writes 

/ it on the drum 

/ jdp clean — this terminates the writing of text onto the drum, 

/ clears buffers, etc. 

/ jdp a set — this initializes things for appending to an 

/ existing text file. The ac contains the field 

/ number where the text starts. Write using 

/ jdp writec and terminate with jdp clean. 



/for reading from the arumf 

/ jdp rset — this initializes things for reading from the drum, 

/ . The field number to be read from is in the ac. 

/ jdp r200 — read one character from the drum into the ac 

/ jdp r400 — read one character from the' drum into the io. 

/ A 14 indicates the end of tex- 



/All routines may be used in nam, lam, or dam, and they 
/save the address mode. - ' „ . -I *,.*, , - • - - y> ? -' * 

/for any reading or writing, a 13 is the end of page character, 



/see D. Thiel for additional information, problems, bugs, etc. 

alio 



/ D. Thiel V^T/JO 

/modified C. landau 22/10/70 
/modified by D. Thiel 1/4/71 
/modified by D. Thiel 18/6/71 



7.1 

ET reading and writing routines 



/call with jdp rset, AC contains field number of text 



/rbl must be a submultiple of 4 00 and a multiple of 40 
rbl=400 . . 

' ? ' *■•■' -1 ' / . '7 /" 



dimension rbf (rbl) 



rset, 

dap rsr 



dac ac 

dio io 

.rar 6s 

add (400-rbl 

dac ddp 

dzm rpc 

lac (isp rpc 

dac rfg /flag for end of file 

law rbf+rbl-1 /last word of buffer 



dap inp 
law rbf 
mta 

lio (340 
law 40 
rsr, ivk . 

hit 
lac ac 
rar 6s" 
add rbf 32 
sub (20000 
dac dne 
lio io 
lac ac 
jmp i rset 



8.1 
ET reading and writing routines 

/call with jdp readch, returns character in io in ivk format 

r2 00, /read into ac 

dac ac 
dio io 
lia 

jdp readch 
lio io 
jmp i r200 

r400, /read into io 

dac ac 
dio io 
jdp readch 
lac ac 
jmp i r400 



rea dch , 

rfg, isp rpc /becomes jmp rp2 after end of file 

jmp inp 
v law i 3 

dac rpc 

idx inp 

dap inq 

sas (lio rbf+rbl 

jmp 9 k 

law rbl 

adm ddp 

lia 

ral 6s 

dap rdr 

law rbf 

dap inp 

dap inq 

mta 

law rbl 
rdr, ivk 

hit 

ZAP 
9k, sub (lio rbf 

add ddp 

sad dne 

jmp rp3 
inp, lio 

cla-~- — y 
Crcl 6s ' 
inq, dio 

sas (77 

sad (13 

dzm rpc 
rpl , lia 

jmp i readch 

rp3, lac (jmp rp£ 

dac rfg 
rp2, law 14 

jmp rpl 



9.1 

ET reading and writing routines 

/et buffer write 

/of 1 must be a submultiple of 400 and a multiple of 40 

/call with jdp wset, ac contains field number of text to be written 

of 1=400 .'> ■ ' ? \ t, , a 

dimension ofb(ofl) .« <'*><TrtS ~ -.*.»* + <* j>y £ 

wset., A ,■->) /* f* 5 i-pfpsj/} 

dac ac '" """ ^ 

dap wrdr 
mta 300 

bpt_ 
rr" bs^- 
-wfld 
add (400-ofl 

dac_ odp_ 

(law i 3. (:: — 
d^c'Cp 
law ofb 
dap ofp 
dap ofq 

(atv77 lac ac 

j ac jtofl cist /last character 

jmp i wset 



10.1 



ET reading and writing routines 



wlOO, 



w300, 



writec, 



wo5 , 



wrd, 

of p-1 , 
ofp, 



ofq, 



v. 



wox, 





dac ac 

dio 95P 
jdp write c 
lac ac 

lio 95P 
jmp i wlOO 



dac 95P 

dio ac 

lai 

jdp writec 

lac 95P 

lio ac 

jmp i w300 



and (77 

dac 96p 

sas (13 

Jmp . 4 

law 77 

sas cist 

jmp force 

3a c (lio ofb+ofl 

sas ofp 

jmp of p-1 

law ofl 

adm odp 

lia 

and (770000 

sas odp 

jmp wo5 

raT~6"s^ 

mta 300 

bpt 
jmp_^+2 

<£al 6s '^ 

d^Twrd 

law ofb 

dap ofp 

dap ofq 

mta 

law of 1 20 

ivk 

hit 

lac 96p 

lio 

rar 6 s 

rcl 6 s 

dio 

isp cp 

jmp ou3 
•• v law i 3 : ~^ -■ - 

dac cp 

idx ofp 

dap ofq 

lac 96p 



/write from ac 



/make sure this is preceded by 77 



/v/rite out buffer 



10.2 

ET reading and writing routines 

dac cist /character just written 

jmp i writec 

ou3, lac 96p 

sas (77 
sad (13 

jmp ofp 
jmp wox 

force, lac writec 

dac acw 
law 77 
jdp writec 
lac acw 
dac writec 

jmp writec 1 



ET reading and writing routines 

Hj — 
/call with jdp clean, finishes up output 3 writes out bufer 

clean, 

dac acl 

dio io 

law ±3 

sas cist /force 13 to precede end of text 

jdp wioo **Wt^*» Hr/iiec. 

lac of p 

add odp 

sub wf Id 

add (20000+of l-[lio ofb] 

dac b77 
law ofb+ofl 
dap ofp 

jdp writec /write out last buffer 

law b77-32 
mta 

lio (340 
law 60 
wrdr, ivk • 

hit 

lac acl 
lio io 
jmp i clean 

b77, 

20400 



12.1 

ET reading and writing routines 

/call with jdp aset, ac contains field number of text to be appended to 
aset, 



ard, 



ardl, 








dac 


ac 


dio 


io 


dap 


ard 


dap wrdr 


, rar ' 


6 s-™ 


dac 


wf Id 


law 


ofb 


mta 




11 o 


1340 


law 


40 


ivk 


• 


hit 




3ac 


wfld 


sub 


(20000 


add 


ofb 32 


dac 


acl 


and 


$[ofl 


lia 




sub 


(ofl 


dac 


odp 


lai 




,ral 


6s-' 


dap 


ardl 


law 


ofb 


mta 




add 


acl 


AMIA 


dap 


ofp 


dap 


ofq 


law 


1 3 ' 


dac 


cp 


law 


ofl 


ivk 




hit 




law 


13 


dac 


cist 


lac 


ac 


lio 


io 


jmp 


i aset 



-1 j ^7/ ?/? 7 



/end of page was last character 



13.1 

ET reading and writing routines 

constants 
variables 
start 



Punch parity routines 
5 February 1972 



Punch parity routines 

/parity generator 

/uses all three registers 

/high 12 bits of register must be clear 

TIAX /or TA3X or TXIA 

sar Is 

sza 

jmp .-3 
SII 
sil 7s 

XVII /or XVIIA etc. 

/garbage left in high 10 bits 



1.1 






.?-« 



«. 



j 









2.1 

Punch parity routines 

/parity generation subroutine for low order 6 io bits 

Ac saved, only low 8 bits of io are significant rest is garbage 

/entry jda pty 

Pty, 

dap ytp /return address 

]aw 2325 /parity bits for io and io bit for ril inst 

rcr 277 

ral 7s Ac contains 0022xx, io 5 2 gggg 

and .-2 

dap . 1 /address patrt includes io bit 

ril .-. /rotates io one plus parity of char 

rar 7s 

rcl 7s /rotate char with parity back : 

lac pty /restore ac 

ytp, jmp . /return 

start 



3.1 

Punch parity routines 

/parity generation subroutine for low order 7 lo bits 

Ac saved, io will contain nothing but char with parity in low 8 bits 

/entry jda pty 

Pty, 

dap ytp /return address 

law 2325 /parity bits for io and io bit for ril inst 

rcr 7s 

ral 7s 

dap . 1 

ril •-• /rotates correct parity bit inl 

rar 7s 

rcr 1 /prepend parity to char 

cli /flush io garbage 

rcl 8s 

lac pty /restore ac 

ytp, jmp . /return 

start 



4.1 

Punch parity routines 

/super short parity generation subroutine for low 7 io bits 
Ac not saved, io bitsO-9 will be garbage 
/entry jsp pty 

/return address 
/parity for io and io bit for ril inst 



/rotates correct parity bit in J 

/rotate char with parity back : 
/return 



Pty, 


dap ytp 




tew 2325 




rcr 7s 




ral 7s 




dap . 1 




n jl # *■• • 




rar 7s 




rcl 7s 


ytp, 


jmp . 


start 





Light pen following routines 
5 February I972 



1.1 

Light pen following routines 

/pen follow - rst - aug 66 

/ts~ only 

/entry? jdp tpn, returns x and y in ac and io resp 

tpn, 

lac x 

lio y 

stf 3 

jmp tpl 
tpl, add inl 

swp 3 

add in2 

swp 
tpl, dpy-i 300 /insert nop before here for non-ts 

szf 3 

jmp tpl 
tpd, sub x 

sar 1 

a dm tx 

lai 

sub y 

sar 1 

a dm ty 

lac inl 

lio in2 

cma 3 

dac in2 

dio inl 

spiVspa i 

jmp tpn 1 

lac tx 

lio ty 

dac x 

dio y 

dpy 300 

szf 3 

jmp i tpn 
fpl, dpy-i 300 

dac tx 

rar 1 

xor (110371 

add (H037I 

lia 

rir 9s 

szf i 3 

jmp fpl 
fpd, lac tx 

lia 

rir 9s 

dac x 

dio y 

dio ty 

jmp tpl 

inl, -1000 

in2, -0 

variables 

constant 

start 



2.1 

Light pen following routines 

/pen follow 

/jdp tpn 

/returns here if can 1 1 f ind pen 

/skips if pen found, coordinates in x,y and A, I 

tpn, 

dzm tx 
dzm ty 
dzm tc2 

lac tsh 
tpl, and (777 

sal Is 

cma 

dac tcl 
tp2, lac ran 

rar 7s 

xor (311071 

add (311071 

dac ran 

scr 9s 

sir 9s 
tsh, rcl 9s 

add y 

swp 3 

add x 

dpy 300 

szf i 3 

jmp tp3 

sar 4s /pen seen 

a dm tx 

lai 

sar 4s 

adra ty 

idx tc2 

law i 17 

ior tcl /don ! t look at more than 17 more points 

dac tcl 
tp3, isp tcl 

jmp tp2 

lac ty 

mul (1 

div tc2 /no. of points under pen 

jmp tp4 /none 

sal 4s 

dac y 

idx tpn /skip return 

law 3007 
dap tsh 
lac tx 
mul (1 
div tc2 
hit 
sal 4s 

dac x 
lio y 

jmp i tpn 

tp4, lac tsh /increase window size 

SAA 
and (777 



2.2 

Light pen following routines 





sza i 






jmp i tpn 


/not found 




adm tsh 






jmp tpl 


/try again 


ran, 


123456 




tx, 







ty, 







x, 







y> 







tcl, 







tc2, 








Random number programs 
5 February ±972 



1.1 

Random number programs 

/simple random number generator Q~ *> «*i*4c > ±* \ J \ 

/good for most applications <*vh*vc fa&% 4o be. *ncr*i*ft*l(Kv / 

ran, 

lac random 

rar 7 s 

xor (311071 

add (3H07I /use carry function 

dac random 

jmp I ran . 



2.1 

Random number programs 

/new improved random number generator 

ran, 

lac rl 
xor (311071 
add (311071 
rar 9s 
dac rl 
lac r2 
xor (355671 
add (355671 
rar js 
dac r2 

xor rl 

jmp i ran /return with number in ac 



Random number programs 



3.1 



/gauss 

/Routine to generate a random number with an approximate gauss ian density. 
/Entry is jdp ran. The density has zero mean and unit variance, with 
/the decimal point assumed between bits 3 and 4. 



ran, 



ml, 



lrn, 



m2, 



law 


1 8 


dac 


set 


dzm 


ras 


law 


i 5 


dac 


ret 


lac 


rda 


rar 


Is 


xor 


(311071 


add 


(311071 


dac 


rda 


and 


(377 


add 


(tbl 


dac 


ads 


lac 


rdm 


ral 3s 


lia 




and 


(177 


spi 




cma 




dac 


inc 


dzm 


rdm 


lac 


rdm 


ral 


3s 


add 


i ads 


dac 


rdm 


idx 


ret 


sma 




jmp 


rnH 


lac 


ads 


add 


inc 


and 


(377 


add 


(tbl 


dac 


ads 


jmp 


Irn 


lio 


ras 


rir 


'ds 


lac 


rdm 


spi 




cma 




dac 


rdm 


a dm 


ras 


isp 


set 


jmp 


rnl 


lac 


ras 


mul 


(122342 


jmp 


i ran 



/Prepare to add 8 random numbers 



/Initialization to form 1 random 15 bit number. 

/Generate a pseudo-random address for the first e 
/to the random number table. 



/Generate a random increment between table entrie 



/Generate a 15 bit random number from 5-3 bit 
/table entries. 



/Complement the number at random. 



/Add to the sum of the previous random numbers 



/unit variance 



Random number programs 
/foandom Number Table 



4.1 



tbl, 




5 


3 


1 


1 


6 


1 


7 


5 


4 








5 


3 


5 


2 


1 


5 


5 


4 


7 


5 


1 


6 


4 


4 





4 


3 


7 


3 


3 


7 


6 


6 


3 


2 


4 


5 


1 


1 


1 


3 


3 


1 


7 


3 





7 


5 


7 


3 


1 


3 


1 





6 


4 


5 


1 


6 


2 





3 





4 


7 


3 


4 


2 


4 


6 


6 





3 


2 


2 





6 


1 


7 


3 


1 


l 


3 


3 





4 


7 


6 


6 


1 


3 





7 


5 





5 


5 













4 

5 
2 

5 
2 

7 

1 

l 

7 

o 

i 

2 

4 

3 
1 

7 

o 



l 

5 
5 
5 



7 

1 
1 

3 
2 

5 
1 
2 
1 
3 
3 
4 
6 
1 
4 
4 

2 

1 

6 
3 

1 

5 

3 
4 

3 
7 
3 



7 

o 

2 

3 

6 
3 
3 
6 
6 
4 
2 
2 

5 
2 



4 
6 
6 
6 
6 
3 
2 
2 
2 

6 

2 

5 

7 
1 

6 

6 

6 

4 

2 

1 

3 
3 

1 

7 
7 

5 
6 
6 
3 
2 
1 
4 
6 



1 

3 
1 
6 
4 
3 
3 

2 

4 

2 
1 
2 
6 

7 

3 

7 

2 

2 

5 
4 

7 

o 

4 

3 

1 

2 

4 
1 
6 

5 
4 

2 

3 
7 
1 

7 

3 
7 
3 

7 

3 
4 

3 
3 
5 
6 
o 

2 
3 



Random number programs 



5.1 



/test random number generator by distribution of sum of successive 
/random numbers. 



go. 



foo, 



bar, 
sea, 



dis, 



n, 

tab, 
tab k/ 



k=1000 

iam 

Ixr (-k 

dzm i tab k 1 

771622 

jmp .-2 

law i 1 

add n 
lio (add 

I 

. 4 

1 

1 

i inir | 



sza 

jmp 
sil 



sar 
jmp 
lai 

ior (sar 
dac sea 
lac n 
cma __ 
dac cint 
dzm tot 
jdp ran 
sar 

adm tot 
Isp cnt 
jmp bar 
lac tot 
sar 9s 
and (777 
add (400 
and (777 
774020 
idx I tab 
szs i 10 
jmp foo 
lac £add 
dac x 
771020 
lac i tab 
ral 9s 
add (add 
lia 

lac x 
lot 307 
771 620 
law 1000 
adm x 
sas (add 1 
jmp dis 
jmp foo 

1 




/determines vertical scale 



/number of terms in s 



um 



ran 







5.2 

Random number programs 

/insert random number generator here, return with 
jmp i ran 

variab 
consta 
start go 



6.1 
Random number programs 

/Test random number generator by correlation of a number with its 
/n-th successor. Put upper and lower limits in TW 0-8 and TW 9-17 
/respectively with signs in TW and 9* Program plots distribution 
/of those random numbers which were preceded on the n-th call before 
/by a random number in the set range. 

go, 



foo, 



bar, 



dis, 



k=1000 


iam 




Ixr 


(-k 


dzm 


i tab k 1 


77I622 


jmp 


.-2 


latVcli 


rcr 


9s 


ral 


9s 


sar 


1 


sir 


1 


dac 


ulm 


dio 


Tim 


771020 


jdp 


ran 


dac 


i buf 


771620 


lac 


n 


775402 


jmp 


..-5 


dzm 


bfp 


Jdp 


ran 


dac 


new 


Ixr 


bfp 


lac 


i buf 


sap 


1 


lio 


11m 


772 613 


jmp 


bar 


lio 


ulm 


772607 


jmp 


bar 


lac 


new 


sar 


9s 


and 


(777 


add 


(4oo 


and 


777 


774 C20 


idx 


i tab 


lac 


new 


Ixv 


bfp 


dac 


i buf 


idx 


bfp 


lio 


n 


773415 


dzm 


bfp 


szs 


i 10 


jmp 


f 00 


lac 


£add 


dac 


X 


771020 


lac 


i tab 



Random number programs 



6.2 



tab, 
ran, 



ral 9s 

add (add 

lia 

lac x 

lot 307 

771620 

law 1000 

adm x 

sas (add 1 

jmp dis 

jmp f 00 

1 

tab k/ 



lac rl 

xor (3HO7I 

add (3IIO7I 

rar 9s 

dac rl 

lac r2 

xor (355671 
add (355671 
rar 7s 
dac r2 
xor rl 
jmp i ran 



/correlation number 



/example of good random number generator 



variab 
consta 
buf , 
start go 



Plotter routines 
5 February I972 



1.1 

Plotter routines 

/move the plotter 

/f and i have signed x and y increments 

/jdp move nam. lam, or dam xr saved 

/does all the 45 degree moves first, then the remaining or 90 degree moves 

/this is a minimum time path 

/can be used to draw perfect 0, 45* or 90 deg lines 

/others will have a bend 

move, 








stf 


5 


taa; 


> 


cma 


5 


dac 


xc 


TIIA<M 


cma 




dac 


yc 


lac 


xc 


TAAP| 


jmp 


• 1 r 


law 


10 


szf 


I 


law 


TII= 




SAA 




TI< 
SAT 




cks 




ril 


7s 


TI>P 


jmp 


.-3 


TAI 




lot 


1111 


isp 


xc 


jmp 


.+3 


law 


i 14 


AMI 


isp 


yc 


Jmp 


.+3 


law 


i 3 


AAIAP 


jmp 


.-14. 


jmp 


i move 
















xc , /x c ount 

yc, /y count 



2.1 
Plotter routines 



dpy plot 

/ss2 up to plot, down to display 

/jdp ini to initialize buffers 

/jdp dsg to gronk dpy 1 s 

/jdp pit to plot a point not originally a dpy 



6000/ 



ini, 









law 


tax 




dap 


ppx 




dap 


• 2 




dzm 


tae 




dzm 






idx 


.-1 




sas 


.-3 




jmp 


.-3 




law 


tay 




dap 


spy 




dzm 


W 




dzm 


py 




jmp 


i ini 


b=100 






tax, 


tax 


b/ 


tay, 


tay 


b-1/ 


tae, 







Pit, 







dsp, 


iot 307 




skp 


i 20 




jmp 


i pit 




dac 


sac 




dio 


sio 


siz , 


sar 


9s 




sir 


9s 


PPX, 


dac 




ppy, 


dio 






lac 


(add-3 




dac 


ds 




law 


tax 




dap 


ppx 




law 


tay 




dap 


ppy 


nex, 


lac 


i ppx 




sub 


W 




dac 


dx 




spa 






cma 






lia 






lac 


i ppy 




sub 


py 




dac 


dy 




spa 






cma 






772610 




swp 






sub 


ds 




sma 





0$5>o l~£7E 



A-I,skp on >0 



2.2 
Plotter routines 

jmp plz 
adm ds 
lac dx 
dac Fix 
lac dy 
dac fly 
lac ppx 
dac spx 
lac £py 
dac spy 

plz , idx ppx 

idx ppy 

sas (dio tae 1 
jmp nex 
lac spx 
dac ppx 
lac spy 
dac ppy 
lac 1 ppx 
dac px 
lac 1 ppy 
dac py 
law 4 
sub ds 
spa 

jdp 1ft 
jdp In 
jdp drp 
dzm f iy 
dzm fix 
jdp In 
lac sac 
llo sio 
jmp i pit 



/dpy gronker 






dsg, . -\ i f *, - 

law \ I* ;^ -W' ,.-jtV 

dap . 1 . Nn" VKm f" 

dii, tec /^ H U . I 

and (760077 -r ^' \v A 

sad (lot 7 y v ^ 

jmp dl2 ni^^ 

di3;, idx dU V** 

sas (lac i 

jmp dU 

jmp i dsg 
d!2, lac dii 

sas (lac dsp 

sad (lac (760077 

jmp di3 

lac (jdp pit 

dac i dil 

jmp di3 



3.1 
Plotter routines 



1ft, 



drp, 



pnd, 



pis, 








lio 


(40 


cla 




dap 


plw 


jdp pll 


jmp 


i 1ft 







law 


600 


dap 


plw 


jmp 


i drp 


dlo 


svl 


lio 


(20 


lot 


1111 


law 


i 14 00 


771312 


jmp 


• -1 


lio 


svl 


cla 




dap 


plw 


jmp 


plw 2 







repeat 2,jdp pll 


jmp 


i pis 



/A+l-*A,skp on 2° 



4.1 



Plotter routines 



In, 

lac fix 
ior fly 
sza 1 
jmp lm 
jdp lin 
jdp pis 
jdp txi 
jmp i In 

lin, 

lio fix 
law 1 
spi 
law 2 
dac lgd 
spi 
cmi 

dio fix 
lio fly 
law 10 
spi 
law 4 
dac smd 
spi 
cmi 
lai 

dio f iy 
sub fix 
spa 

jmp 111 
lac fix 
dac f iy 
dio fix 
lac smd 
lio lgd 
dac lgd 
dio smd 

111, lac fix 

cma 

dac npt 
cma 
sar 1 
dac nt 

tis, lac fiy 

adm nt 
sub fix 
cli 60 
spi 

jmp . 3 
dio nt 
lac smd 
ior lgd 
swp 
jmp i lin 

txi, 

isp npt 



4.2 

Plotter routines 

jmp tis 

jmp i txi 

lm, lio (1 
jdp pll 
lio (2 

jdp pll 
jmp i In 



Plotter routines 



pll, 

plw, skp i 



constants 
variables 
start 



/jmp .-1 



jmp i pH 



5.1 



jmp pnd 

lot 1111 

law i 240 

/771312 /A+l-^,skp on 2° 



add (1 
spa 
jmp .-2 



6.1 



Plotter routines 



/lineplot, 11 sept ±969 

/"getccm" will assign the plotter 

/It also helps to turn on the Calcomp, etc. 

/to start, line up your origin, then dzm px and dzm py 

/the macros M drawto x ,y* will move the pen from its 

/present position to (x,y), where x and y are addresses of floating 

/point numbers. The x and y scales given here are for 

/1.0 at seven inches. You can change them during your 

/run. Watch out for going off the graph. 



/"penup* and "'pendwn 1 * do the obvious thing 



def ine 



sea lex, 
sea ley, 



draw to x,y 

load y 

mulby sea ley 

jdp fix 



dac f bl 

load x 

mulby sea lex 

jdp fix 



lio f bl 

jdp pin 

terminate 

205657 
205657 



n 



$50 L&T££ 



/1400.0 scales 1.0 as seven inches 
/other scales are possible 



penup=jdp pup 
pendwn=jsp pdn 



define 



getccm 
law f lexo 
arq 
bpt 
terminate 



q2 



7.1 

Plotter routines 



/plot from (px,py) to (AC,IO), then update (px,py) 

/everything saved 

/runs in nam, iam, or dam 

/coordinates are signed integer, step size is .005 inches 

pin, /jdp 

dac tac 
dio tio 
sub px 
dac dx 
adm px 
lai 

sub py 
dac dy 
adm py 
lac dy 
lio dx 
sil 2s 
scl 2s 
law 12 
A~IA 

and (±7 
dac dag 
lac dx 
lio dy 
spa 
cma 
spi 
cmi 
dac dx 
AMI< 
jmp .44 
dio dx 
lio (-14 
A-»IA | 
law i 3 
and dag 
dac str 
dio dy 
dzm sds 
dzm Ids 
plb, lac Ids 

sad dx 
jmp ret 
idx Ids 
mul dy 
div dx 
jmp ret 
lio str 
sas sds 
lio dag 
dac sds 



8.1 
Plotter routines 



ret, 

pup, 6 ^ /pen up jdp 

pc2, 

pdn, dac pup /pen up jsp 

Pit, 



jdp 


pit 


jmp 


plb 


lac 


tac 


lio 


tio 


jmp 


i pin 







lio 


(4o 


jdp 


pit 


jmp 


i pup 


dac 


pup 


lio 


(20 


jmp pc2 







3ai 




cks 




ril 


7s 


spi 


i 


jmp 


• -3 


lia 




lot 


llll 


jmp 


i pit 





















































tac, 
tio, 
dag, 

str , 

Ids, 

sds, 

dx, 

dy, 

px, /where pen is now 



9.1 



Plotter routines 



/Incremental line plotting routine 
/call with jdp In 
/dx in A, dy In I 



In, 



lne, 



lnf , 



Pl» 





A-*3X 
law 4 
X-»DC>P 
CXX | 
law 10 
dac lna 
TII>P 
CII| 
SAA 
SAA 

dac lnb 
XMI< 
jmp lne 
X-KDC 
xor lna 
dac lna 
CXXA 
dac lnd 
dio lnc 
sar Is 
SXX<= 
jmp i In 
sub lnc 
sub lnd 
lio lnb 
spq 

jmp .+3 
add lnd 
lio lna 
dac lnf a 
lai 

jdp pi 
lac lnf a 
jmp lnf 




jmp .-3 
lia 

lot 1111 
jmp i pi 



/CAX for Howell standard 



Qfr'pOl.e'Xe 




10.1 
Plotter routines 

/absolute line plotting routine 

/call with jdp In 

/x in A , y in I 

/absolute pen positions stored in px, py 

m, o nb^oLB^E 

sub px [y is ^ 

<?et pert"* 
library 



lne, 



lnf , 



Pi, 



pi, 








sub 


px 


swp 




sub 


py 


TAXP 


i 


adm 


py 


TIIA 


adm 
law 


F 


X-»IX>P 


CXX " 




law 


10 


dac 


lna 


TII>P 


CII| 


1 


SAA 




SAA 




dac 


lnb 


XMI< 


jmp 


lne 


X->IX 




xor 


lna 


dac 


lna 


CXXA 




dac 


lnd 


dio 


lnc 


sar 


Is 


SXX<= 


jmp 


i In 


sub 


lnc 


sub 


lnd 


lio 


lnb 


spq 




jmp 


.+3 


add 


lnd 


lio 


lna 


dac 


lnfa 


3a 1 




jdp 


Pi 


lac 


lnfa 


jmp 


lnf 







cks 




ril 


7s 


jmp 


.-3 


lia 




lot 


1111 







cks 




ril 


7s 


spi 


i 


jmp 


.-3 


lia 





10.2 
Plotter routines 

lot 1111 

jmp 1 pi 



11.1 

Plotter routines 

/character plotting routines on eight pages 

/txp 

/text plotter 

/define ori and siz, then bring pen to location 

/jdp txp, followed by the text, ending with a • 



txp, 



tpx, 11 o i txp 

idx txp 






rpf 
dio 47t 

nam ^ *- .. > ,■--,- ;, n 



/ 
€ C 



lac (607600 | ^py 

rcl 6s / 

sa s ( la i 

jmp . 6 

lio 47t 

lpf 

lac (nop 

dac chl-1 

jmp i txp 

sad (swp 

jmp tpx 

dac 48t 

dio 48s 

jdp cpl 

3a c (jmp ch2 

dac chl-1 /so get right place for next lette: 

lac 48t 

lio 48s 

jmp tpx 2 



12.1 

Plotter routines 

/character plotting subroutine 

/before start, need to assign calcomp, law flexo q2, arq, hit 

cpl, .-• /use jdp cpl with the character in A (in concise 

/code) and the orientation in ! ori 1 and the size code in ! siz f 

/size code is -1 for eighth inch, -2 for quarter inch, 

/up to -20 (-16.) for two inch (max) 

/txp changes to jmp ch2 

come in otherwise, need to restore 

/communicate with pin 





dac 


ch 




nop 




chl, 


rpf 






dio 


svf 


/runs 


in nam, t 
nam 


30 if 




lac 


px 




dac 


xO 




dac 


x2 




lac 


py 




dac 


yo 




dac 


y2 


ch2, 


law 


77 




and 


ch 




add 


(dsp 




dac 


get 




lio 


i get 


cas, 


skp 


i 600 




ril 


9s 




cla 






rcl 


9s 




sub 


(6 




spa 






jmp 


ttbl 




add 




dac 


get 




lio 


i get 




cla 






rcl 


4s 




add x2 




dac 


xl 




cla 






rcl 


5s 




sad 


(37 




jmp 


spo 




add 


y? 




dac 


£i 


nO, 


dzm 


t 




lac 


yi 




sub 


yo 




sza 


i 




jmp 


nxl 




and 


(4000 




sza 






law 


i 4 




add 


(10 




dac 


t 



/special origin 



12.2 
Plotter routines 

nxl, lac xl 

sub xO 

sza i 

jrap nyl 

and (400000 

ral Is 

add (1 

adm t 
nyl, lac t 

sza i 

jmp go 

ior (40 
jda pit 

jmp nO 



13.1 



Plotter routines 

spo, Idx get 

lac i get 
add x2 
dac xl 
idx get 
lac i get 
add y2 
dac yl 
Idx get 
jmp nO 

go, 3ac (add 

adm get 



law 300 
and t 
sza 1 
jmp nx2 
law 4 00 
and t 
rar 6s 
cma 

add (10 
adm pit 

nx2, law 200 

adm t 
and (300 
sza i 
jmp ny2 
law 4 00 
and t 
rar 8s 
add (1 
adm pit 

ny2, law 37 

and t 
cma 
dac t 
jsp plt+1 
isp t 
jmp .-2 
jmp go 



pit, 



lio 1 get 

spa 

ril 9s 
cla 
rcl 9s 

sad (777 

jmp dun 

sad (776 

jmp 1 cpl 

dac t 

and (40 

rar Is 

cma 

add (40 

dac pit /pen up/down 



Plotter routines 



13.2 



pll, 



lrx, 



dap 


plx 


lac 


siz 


dac 


sz2 


llo 


pit 


law 


1 


and 


ori 


sza 




jmp 


lr 


law 


2 


and 


ori 


sza 




jmp 


ud 



/size code, from main 



/orientation, from main prog 



14.1 
Plotter routines 

udx, lot 1111 

/figure actual location for pin 

/test -y "bit 





law 


i 1 




rir 


Is 




spi 






a dm 


py 




law 


i 1 




rir 


2s 




spi 






adm 


PK 




ril 


Is 




spi 






idx 


py 




rir 


2s 




spi 






idx 


px 




lac 


pit 




xor 


p32 




and 


(60 




sza 






law 


i 2000 




sub 


I C~»L~mC~m f 




dac 


P12 




isp 


p32 




jmp 


.-1 




lio 


pit 




dio 


pl2 




isp 


sz2 




jmp 


pll 




rir 


Is 




spi 






idx 


xO 




rir 


Is 




law 


i 1 




spi 






adm 


xO 




rir 


Is 




law 


i 1 




spi 






adm 


yo 




rir 


Is 




spi 






idx 


yo 


plx, 


jmp 


• 


P12, 







lr, 


rcr 


2s 




ral 


Is 




rcl 


Is 




rar 


2s 




rcr 


3s 




ral 


Is 




rcl 


Is 




rar 


2s 




rcl 


Is 




ral 


Is 




rcl 


2s 




jmp 


lrx 



/update position 



14.2 



Plotter 


routines 


ud, 


rcr 4s 




ral 2s 




rcl 2 s 




rar 3s 




rcl Is 




rar 2 s 




rcl Is 




jmp udx 



Plotter routines 



15.1 



gl, 



bks, 



cr 



tab, 



dun, 



add 


[xf 1 


dap 


jmp 


. 


jmp 


dun 


jmp 


tab 


jmp 


cr 


jmp 


bks 


ZAP 




law 


600 


dap 


cas 


jmp 


dun 2 


law 


i 22 


adm 


x2 


jmp 


dun 2 


law 


i 40. 


adm y2 


dzm 


x2 


jmp 


dun 2 


law 


CLcL % 


add 


x2 


mul 


(220. 


div 


hit 




add 


(1 


mul 


(220. 


div 


(l 


hit 




dac 


x2 


law 


22. 


adm 


x2 


lio 


svf 


lpf 




jmp 


i cpl 



/space 

/tab 

/car. ret. 
/backspace 
/upper case test 



/advance to next character 



/exit 



16.1 



Plotter routines 



dsp, 



6 
6 
6 
6 
6 
6 





ql-tbl 6 
[qu2-tbl 
q3-tbl 6 
q4-tbl 6 
q5-tbl 6 
q6-tbl 6 
q7-tbl 6 
q8-tbl 6 
q9-tbl 6 







qO-tbl 6 
[qqq-tbl 
qs-tbl 6 
qt-tbl 6 
qu-tbl 
qv-tbl 
qw-tbl 
qx-tbl 
qy-tbl 
qz-tbl 


[qeq-tbl 





1001 



[qub-tbl 
qj~tbl 6 
qk-tbl 6 
ql-tbl 6 
qra-tbl 
qn-tbl 
qo~tbl 
qp-tbl 
qq-tbl 
qr-tbl 



iqpl-tbl 

'qrb-tbl 
qvb-tbl 
.qlb-tbl 



qa-tbl 

qb~tbl 

qc-tbl 

qd-tbl 

qe-tbl 

qf-tbl 

qg-tbl 

qh-tbl 

ql-tbl 6 

5005 



/space 
6]xl000 q2-tbl 6 



6]xl000 qsl-tbl 6 



6]xl000 qcm-tbl 6 

/tab 
6 ]xl000 qcd-tbl 6 



6 
6 
6 
6 
6 
6 



6]xl000 qmn-tbl 6 

61x1000 qrp-tbl 6 
6]xl000 qob-tbl 6 
6]xl000 qlp-tbl 6 



6 
6 
6 
6 
6 
6 
6 
6 



/lower case 



16.2 
Plotter routines 

[qtm-tbl 6 klOOO qpe-tt>l 6 /period ,x 

4004 /upper case 

3003 /backspace 



2002 /car. ret. 



Plotter routines 



17.1 



tbl, 
/numbers 

103743 

qi, 

2ooo46 

q2, 

740457 

03, 

004744 
544710 

5^0270 
777000 

004744 

014053 
7 44 777 

Q7, , 

542644 

q8, 

214544 

00^744 
053777 



045143 

403270 

241157 

047144 
043777 

406642 

047144 
744644 



057643 
777000 

644744 
344447 

047144 



/upper case numbers 

qu2, /single quote 

351257 041657 
/punctuation 

qsl, /slash 

000244 157245 

qqq, /question mark 

403441 541641 

£46145 

£4 434 5 

qpe, /period 

340041 141241 

qtm, /times sign 

106154 414754 

qcm, /comma 

440442 341241 

qeq, /equals sign 

111054 206454 

qcd, /center dot 

352041 141241 

qub, /under bar 

777000 -6 

qob, /overbar 

777000 -6 

qvb, /vertical bar 

777000 7 

qlp, /left paren 

500346 254146 



262343 
544777 
244344 

uttj r r 

541643 

244344 
544447 

542643 



047144 
260344 



777000 

777000 

741041 

444 54 ^ 

341441 

777000 

141041 

777000 

341441 

-4 

34 

-4 

777000 



445543 

447544 
144244 

541642 

453254 
344260 

542643 



244344 
744047 

447544 



141241 
644777 
541641 

741642 



662777 

642777 
344447 

543057 

057777 
144047 

542643 

144244 
777000 

644744 



341203 



741777 



544777 



541641 


741776 


OOO066 


776000 


000066 


776000 


000277 


776000 



I?6l46 $$& paren 777000 

qlb, /left bracket 

54.0447 270047 777000 



17.2 

Plotter routines 



777000 



qrb, 


/right bracket 


200047 


270447 


qmn, 


/minus sign 


114054 


777000 


qpi, , 


/plus sign 


114054 


306654 



777000 



Plotter routines 



18.1 



/upper case letters 



qa, 
000263 

000^70 
214052 

qc, 

746642 

777000 

qd, 

OOC270 

qe, 

000270 

Of , 
OOO270 

qg> 

510045 

642777 

qh, 

000270 

qi, 
200046 

qj, 

006643 

qk, 

00C270 

CB0670 

qm, 

000270 

qn, 

OOC270 

qo. 

004260 

qp> 
000270 

qq, 

004260 

oo4 o44 

qr, 

000270 

qs, 

004744 

047744 

C80056 

qu, 

030664 

qv, 

C30661 

qw, 

CB0665 

qx, 

000245 

qy> 

340256 

qz, 

740457 

xO, 



145045 

052744 
777000 

544447 



053744 

057514 
057514 

644544 

614057 

403270 

743044 

6I6I56 

057777 
747147 

604757 
144 047 

053744 

144047 
745777 

053744 

047144 
777000 

407670 
744o47 
747147 
743144 
156245 
347243 
245I57 




745663 
644544 

344260 

660544 
403054 
403054 
447344 

214670 
403046 
143265 
514754 

670777 
605270 

744660 
644544 
744660 

644544 
244344 

777000 
144264 

261 777 
744143 
416645 
016643 
244457 



317057 
744644 

144047 

453777 
514057 
777000 
260144 

777000 
777000 
405052 
77700Q 



777000 

265777 
756645 
547777 
777000 



777000 
544452 

744642 



777000 



047744 



777000 



777000 




P rTTT / 


344777 


453777 




544447 


344102 


453003 


754777 


447344 


244144 



777000 



/where pen is now 



yo, o 

x2 , Adhere next letter begins 



18.2 
Plotter routines 

ys, o 

constants 
variables 



Scope routines 
5 February ±97^ 



=4> 



1.01 

/hack test and/or demo program for sdp 
test, e2m 

ckn 
lai 

ckn 100 
LZI 10. 

TIX , 

ckn 200 /' 

1ZI 1 

XVII 

jdp sdp 

jmp test 

/twos mode version of . \ r i « LI* 

/scope decimal print <^r^r ^ U^ */**) a-,W o^ pov>»ic 

///P ^dp „. _ We #» CpJp-f 1 ) vhcUv- fcl****. 

/I0 contains scope coordinates l r ' 

/ x, low 9 bits y, high 9 bits 

/AC number to be displayed Scop€ 

(Aec'i^na 



-ClA/R 
73(%2£ 



size=2 


/power of 2 


sdp, 









dio 


sdz 


sda, 


cliVswp 




div 


• 1 




12 






dac 


s"dy 




swp 






LZA 


1 




add 


(sdt 




dap 


sdc 




dzm 


sdx 




lan 


21 


sdb, 


dac 


sdw 




dac 


sdv 


sdc, 


lio 


• 




idx 


sdc 


sdd, 


lac 


sdx 




sub 


(400001+6007xsize 
(40000CV7007xsize 




and 




spa 






add 


(400000+lOOlxsize 




dac 


sdx 




LRI 


1 




spi 


i 




jmp 


sde 




dio 


sdu 




add 


sdz 




lia 






RRA 


9. 




iot 


207 




lio 


sdu 


sde, 


isp 


sdw 




jmp 


sdd 




lan 


22 




sas 


sdv 




jmp 


sdb 




lan 


6xsize 




add 


sdz 



1.02 





and (776777 








dac sdz 








lac sdy 








sza 








jmp sda 








jmp i sdp 






sdt, 


±75014 


30137 


/o 




27 


360000 


A 




305215 


31143 


/2 




105014 


31133 


/3 




36100 


217704 


A 




137114 


631130 


/5 




175H4 


231131 


/6 




3610 


620501 


/7 




155H4 


231133 


/8 




15114 


231137 


/9 


variab 








consta 








start test 









Scope routines 



1.1 



/scope decimal print 

/jdp sdp 

/IO contains scope coordinates 

/ x, low 9 bits y, high 9 bits 

/AC number to be displayed 



size=l 

sdp, 
sda, 



/power of 2 

_ 

dio sdz / 

cliVswp 



A- 



^+y 



\^ 




/,„ ;-f o<*«-l A'rptp C^JqA 



J.i/11 



sdb, 
sdc, 
sdd, 



add (sdt 
dap sdc 
dzm sdx 
U/\ tew-i 21 
dac sdw 
dac sdv 
lio . 
idx sdc 
lac sdx 



sT 



tfocCso 



sde 



l#A*t 



j.y^ 



sub i400000+6007xsize 

and (40000qgf007kslze v - |K 

spa 

add (400000+lOOixsizefil- juM- 
dac sdx . 

ril 1 " tRj 1 

spi i / ~"~ 

jmp sde 
dio sdu 
add sdz 
lia 

■rar--9s~ R^A ^. ,- 

lot 207 
lio sdu 
isp sdw 
jmp sdd 

sas sdv 

jmp sdb 

law -i 6xsize 

add sdz 

and 1 - 100Q — 1 1 U~7 11 

dac sdz 

lac sdy 

sza 

jmp sda 

jmp i sdp 



-K.. 



+v»»3 



sdt, 



175014 
27 

305215 
105014 

36100 



30137 
360000 

3H43 
31133 
217704 



A 
A 



Scope routines 



variab 
consta 
start 



1.2 



137H4 63H30 A 

175114 231131 /b 

3610 620501 /7 

155114 231133 /8 

15114 231137 /9 



Scope routines 



2.1 



/simple line display 

/jdp dsp with coordinates of endpoints in dxl , dyl , dx2 , dy2 
/coordinates must be right justified, magnitude<1000 
/dxl , dyl , dx2 , dy2 are saved 



dxl, 
dyl, 
dx2, 

dy£, 

dx3, 
dy3, 

dsp, 



/5/-T 



dsi; 



dn, 













lac dx2 

sub dxl 

TAX 

lac dy2 

sub dyl 

dac dy3 

spa 

cma 

15CIX>= 

cmi 
AMI>= 
swp 
sir 1 
A+II 
sir 2s 
dio dn 
TXA 

CUC4 - 

jmp dsl— > 

mul (200 
div dn 
hit 
sal Is 
d ac dx3 

tSc dy^ 

iul_{200 
div dn 
hit 
sal 1 

dac dy3 
lac dxl 
lio dyl 
rcl 8s 
lot 207 

SXXflB ^ 

jmp i dsp 
swp 

add dy3 
swp 

add dx3 
jmp dsl 




/x coordinate of endpoint 
/y coordinate 



/could be modifed to enter with 
/x2 and xl in A I 



/change spacing here 



U^ "^^ 



VVN, 



\ d 



r 



/no, of points -1 



Scope routines 



&*- 









3.1 



/truncating line display with fixed window 

/Coordinates in xl, yl, x2 , y2, right justified, up to 17 bits 
/Coordinates are truncated down to magnitude <^ 777 (size of scope face) 



dsp, 
sqn, 



sqO, 



sq2, 








iara 




law 


tab 


dac 


xrl 


elf 


6 


Ixr 


xrl 


bam 




law 


i 777 


szf 


6 


cma 




lia 




add 


i 3 


swp 




add 


i 4 


szf 


6 


cmaVcmi 


taam| 


cla 




TIIM| 


cli 




AVI<M 


jmp 


lgb 


AAI>P 


jmp 


sq2 


spi 




sxx , 


1 


swp 




dac 


11 


AMIP 


L 


dac 


12 


lac 


i 4 


sub 


i 3 


mul 


11 


div 


12 


hit 




adm 


i 3 


lac 


i 1 


sub 


i 


mul 


11 


div 


12 


hit 




adm 


i 


iam 




szf 


6 


jmp 


.+3 


stf 


6 


jmp 


sqO+1 


law 


3 


adm 


xrl 


sas 


(tab+6 


jmp 


sqO 


lac 


x2 


sub 


xl 


sal 


6s 



/could be nam 



/forget it, can't be seen 
/no adjustment 



3.2 
Scope routines 





lia 




lac y2 




sub yl 




sal bs 




dac lg3 




spa 




craa 




spi 




AMIAX | 




A+IAX 




dac lg9 




lai 




scr 7s 




div lg9 
jmp lg4 






dac lg7 




lac lg3 




scr 7s 




div lg9 


lg3, 







dac lg3 




1KXA 




sar 8s 




CAX 


lg4, 


lio yl 




lac xl 




rcl 8s 


3g5, 


dpy-i 300 




add lg7 




swp 




add lg3 




swp 




SXX>P 




jmp lg5 


lg6, 


rail r\o^m - 
jmp i dsp 


lg7, 





3§9, 
xl, 






yi, 





x2, 





y2, 





tab, 


xl 




yi 




xl 


• +6/ 




constants 




variables 




start 





WO-tj 4^v5» .sUft-JVck W oUWWo^ 



x2 xl 

y2 yl 

x2 xl 



4.1 



Scope routines 



/line drawer with variable window 

/takes beginning coordinates in (xl,yl), end in (x2,y2), these are 
/destroyed. Coordinates are right justified, up to 17 bit magnitude. 
/Coordinates are truncated according to the limits in upr, lwr, 1ft, 
/and rgt, and the result is displayed. Goes to fgt in base mode if 
/no part of line is visible. Displays line once and goes to lg6 if 
/it is visible. Limits should not exceed + 1777 (size of scope face) 



sqn, 

sqO, 

sql, 
sq3, 



sq2, 



11 



lam 

law upr 

dac xrl 

clavclf i 

dap sql 

dap sq3 

lxr xrl 

bam 

lac i 10 

sub 

lia 

lac 

sub 

szf 6 

cmaVcmi 

TAM| 

cla 

TIM J 

cli 

AVI<M 

jmp fgt 

AAI>P 

jmp sq2 

spi 

sxx 

spi i 
swp 

dac 16 
AMIA 
dac 17 
lac i 11 
sub i 10 
mul 16 
div 17 
hit 

adm i 10 
lac i 6 
sub i 5 
mul 16 
div 17 
hit 
adm 
iam 
stf 6 
idx sql 
sas (sub 2 
jmp sq5 
law 3 
adm xrl 
sas (upr +6 
jmp sqO 



i 5 



/could be nam 



/forget it, can't be seen 
/no adjustment 



4.2 



Scope routines 



lg3, 



^4, 

365, 



lg6, 
fst, 

3g7, 

te9, 
xl, 

yi> 

x2, 

upr, 
Iwr, 



lio xl 


sil 


8s 


dio 


xl 


sir 


2s 


lac 


x2 


sal 


6s 


A Mil 


dio 


1S7 


lac 


y2 


sub 


y 1 


sal 


6s 


dac 


lg3 


spa 




craa 




spi 




cmi 




A+IA 


L 


dac 


lg9 


lai 




scr 


7 s 


div 


lg9 


jmp 


lg4 


dac 


lg7 


lac 


lg3 


scr 


7s 


div 


lg9 







dac 


lg3 


lac 


lg9 


sar 


8s 


craa 




dac 


lg9 


lio 


yi 


sil 


8s 


lac 


xl 


dpy- 


-i 300 


add 


1&7 


dac 


xl 


lac 


lg3 


A +11 


isp 


lg9 


jrap 


lg5 


jmp 


4 


jmp 


4 
































1777 


-1777 







Wr 


xl 




x2 




xl 





4.3 
Scope routines 

yi 
y2 

xl 

x2 

xl 
.+6/ 
consta 
variab 
start 



Scope routines 



5.1 



/character display 

/upper and lower case 

/jdp cd, character in ac [12-17 ] 

/ac,io,xr lost 

/nam or iam 

/char 16 =>reset pointer to upper left corner of scope 

/char 14 =>pointer => 

/code word format 

/ L7 Ll4 |R3 RlO R17 

/ L6 1^.3 R2 R9 R16 

/ L5 112 Rl R8 R15 

/ 1/1 Iil RO R7 Rl4 

/ 13 HO H7 R6 R13 

/ 12 L9 H6 R5 R12 

/ Ll L8 145 R4 Rll 



cd, 



cdl, 
cc£ , 

cd3, 



cd4, 



cd5, 





and (77 

sal 1 

add cds 

dap cd4-l 

cla 

dap cd6 2 

law cd5-l 

dap cd4 

xct cd4-l 

spa 

jmp c9d 
and (376000 
TAXA 
lio cdy 
A+XX< 
jmp cd3 
lac cdx 
dpy-i 200 
law 1000 
A +11 

lac (376000 
AAXXA = 
jmp cd2 
law 1000 
adm cdx 
idx cd4 
lac . 
xct . 
jmp cdl 
idx cd4-l 
xct cd4-l 
rcr 4s 
jmp cdl 

ral 7s 
TAI| 
ral 3s 
rar 8s 
jmp cd6 



cd6, 



law 2000 



Scope routines 

adm cdx 
law . 
adm cdy 
jrap i cd 

cd7, ^c (77777 /tab 

ior cdx 
add (7001 
dac cdx 
jmp i cd 



5.2 



Scope routines 

cdc, ZAP /lower case 

/upper case 



6.1 



cds, 
cdx, 

cdy, 

c9d, 



c8d, 



ZAP 




law 


200 


add 


(cd9 


dac 


cds 


jmp 


1 cd 


cd9 




507000 


300000 


sas 


(add 


jmp 


c8d 


idx 


cd4-l 


dap 


. 1 


Jmp 


• 


law 


3000 


dap 


cd6 2 


cma 




a dm 


cdy 
cd4-l 


xct 


jmp 


cdl 


lac 


(-14000 


a dm 


cdy 


lac 


(507000 


dac 


cdx 


jmp 


i cd 


lac 


(300000 


dac 


cdy 


jmp 


cdB+2 



cd8, lac (-14000 /c.r 



ini, lac (300000 /16 - Initialize 



Scope routines 






cd9, 





/space 


27 


740000 


A 


305214 


462306 


/2 


105014 


462266 


/3 


36100 


437610 


/4 


137H4 


462261 


/5 


175114 


462262 


/6 


3610 


441203 


n 


155H4 


462266 


/8 


15H4 


452236 


/9 


add 


jmp 1 cd 




add 


jmp i cd 


/stop 


51253 


305010 


/pointer 


add 


jmp 1 cd 




add 


jmp Ini 


/initialize 


add 


jmp 1 cd 




175014 


60276 


/o 


100200 


401002 


// 


221245 


211000 


/s 


10764 


211000 


A 


171004 


37000 


/u 


704 04 


10034 


/v 


171003 


20074 


/w 


210501 


12104 


A 


417104 


417600 


/y 


211445 


223104 


/z 


add 


jmp I cd 




500200 





/, 


add 


jmp 1 cd 




add 


jmp i cd 





7.1 



add 
add 



jmp cd7 

jmp i cd 



/tab 



Scope routines 



8.1 





501004 

3762 02 

1774 

370047 

370040 

161044 

776IIO 

414110 

370040 

add 
add 
20100 
1012 

2010 
342 

add 

161044 

377104 

161044 

161044 

16124 5 

21760 

415114 

37 61 00 

7 

add 

4 

add 

add 

add 

add 



400000 


/• 


217200 


/J 


421000 


A 





/l 


401170 


/m 


236000 


/n 


216000 


/o 


441400 


/P 


477640 


/q 


202000 


/r 


jmp 1 cd 




jmp i cd 




402010 


/- 


107000 


/) 


40201 


/'" 


120200 


/( 


jmp 1 cd 




237100 


/a 


4l4ooo 


A> 


212000 


/c 


237600 

226000 


/d 


/e 


44o4oo 


/f 


457600 


/g 


434000 


A 


500000 


/i 


jmp cdc 







/. 


jmp cdc 1 




jmp 1 cd 




jmp i cd 





jmp cd8 



9.1 



Scope routines 












/space 


30 


600 


/ M 





140000 


/ f 


4olo 


4o4ol 


r 


104422 


110434 


/O 


70402 


10034 


A/ 


70020 


100434 


/A 


20242 


120200 


/< 


1012 


105010 


/> 


10027 


74o4o4 


/t 


add 


jmp i cd 




add 


jmp 1 cd 


/stop 


51253 


305010 


/pointer 


add 


jmp i cd 




add 


jmp Inl 


/initialize 


add 


jmp i cd 




20102 


507010 


/■> 


4015 


42206 


/? 


115114 


462262 


/s 


2017 


740201 


/T 


177004 


20077 


/u 


76404 


10037 


At 


177003 


20077 
405I43 


/W 


306240 


A 


6047 


401003 


ft 


303214 


461303 


/z 



add 

50241 

add 

add 

add 

add 



jmp i cd 
205024 
jmp i cd 
jmp i cd 
jmp cd7 
jmp i cd 



/= 



/tab 



10.1 



Scope routines 






6oloo4 


20100 


/ 


101004 


20077 


/T 


376101 


210501 


/K 


377004 


20100 


/L 


376020 


200577 


/M 


376040 


404177 


/N 


175014 


60276 


/o 


376110 
175015 


442206 


/P 


50336 


/Q 


376111 


452306 


M 


add 


jmp 1 cd 




add 


jmp i cd 




20103 


702010 


/+ 


1014 


77600 


/] 


7 


740000 


/ 


1774 


60200 


/( 


add 


jmp 1 cd 




370221 


44574 


/A 


377114 


462266 


/B 


175014 


60242 


/c 


377014 


60276 


/D 


377H4 


462301 


/E 


376110 


442201 


/F 


175015 


64262 


/G 


376100 


402177 


/H 


1017 


760200 


/I 


add 


jmp cdc 




104240 


405042 


/x 


add 


jmp cdc 1 




add 


jmp 1 cd 




add 


jmp i cd 




add 


jmp cd8 





11.1 

Scope routines 



variao 
consta 
start 



Scope routines 



12.1 



/character display 

/upper and lower case 

/jdp cd, character in ac[12-17] 

/ac,io,xr lost 

/nam or iam 

/char 16 =>reset pointer to upper left corner of scope 

/char 14 =>po inter => 

/code word format 

/ L? 114 |R3 RIO R17 

/ L6 H3 R2 R9 R16 

/ L5 H2 Rl R8 Rl5 

/ L4 111 ]R0 R7 Rl4 

/ 13 L10 117 IR6 R13 

/ 12 L9 116 R5 R12 

/ II L8 115 R4 Rll 
.ds=lOOO 



/instructions become lac ( f s. 

define amacro w 

repeat ifp 7777 -w, law w 
repeat ifm 7777~w,lac (w 

termin 



/character size control. If size>llll octal, some law 



cd, 



cdl, 
cd3, 



cd4, 





and (77 
sal 1 
add cds 
dap cd4-l 
law 

dap cd6 2 
law cd5-l 
dap cd4 
xct cd4-l 
spa 

jmp c9d 
and (376000 
TAXA 
lio cdy 
A-£CX< 
jmp cd3 
lac cdx 
lot 207 
amacro . ds 
A +11 

lac (376000 
AAXXA = 
jmp cd2 
amacro .ds 
adm cdx 
idx cd4 

1B.Q, . 

xct . 
jmp cdl 
idx cd4-l 
xct cd4-l 
rcr 4 s 
jmp cdl 



JL.CZL (—. 

Scope routines 
cd5, 



cd6, 



ral 


7s 


tai; 


1 


ral 


3s 


rar 


8s 


jmp 


cd6 


amacro .dsx2 


adm 


cdx 


law 


• 


adm 


cdy 


jmp 


i cd 


lac 


(77777 


ior 


cdx 


add 


(7001 


dac 


cdx 


jmp 


i cd 



cd7, tec (77777 /tab 



Scope routines 



13.1 



cdc , 


ZAP 




law 200 




add (cd9 




dac cds 




jmp i cd 


cds, 


cd9 


cdx, 


507000 


cdy, 


300000 


c9d, 


sas (add 




jmp c8d 




idx cd4-l 




dap . 1 




jmp . 


c8d, 


/ amacro .dsx3 
\ dap cd6 2 






v cma 




adm cdy 




xct cd4-l 




jmp cdl 


cd8, 


lac l~.dsxl4 




adm cdy 


cd8 2 , 


lac (507000 




dac cdx 




jmp 1 cd 


Inly 


lac 1300000 




dac cdy 




jmp cdo-f2 


old, 


amacro . dsx7 




adm cdx 




jmp i cd 



c6d, 



/lower case 
/upper case 



/push bottom of 5x7 matrix down to 
/display symbols ,ypqgj 



/ \ i A 



\ 



/c.r. ^ 

/c.r. with no line feed [Jb) 

/lb - initialize 
/space 



repeat ifp 7777-«dsx7> law i .dsx7 
repeat ifm 7777-. dsx7> lac (-.dsx7 
adm cdx 

jmp i cd 



/backspace 



Scope routines 



14.1 



cd9, 



add 

27 

305214 

105014 

36100 

137H4 

175114 

3610 

155114 

15114 

add 

add 

51253 

add 

add 

add 

175014 

100200 

221245 

10764 

171004 

704 04 

171003 

21050I 

417104 

211445 

add 

500200 

add 

add 

add 

add 



jmp c7d 

740000 

462306 

462266 

437610 

462261 

462262 

441203 

462266 

452236 

jmp i cd 

jmp 1 cd 

305010 

jmp i cd 

jmp ini 

jmp i cd 

60276 

401002 

211000 

211000 

37000 

10034 

20074 

12104 

417600 

223104 

jmp i cd 



jmp i cd 

jmp i cd 

jmp cd7 

jmp 1 cd 



/space 

/I 

/2 

/3 

/4 

I 

/stop 

/pointer 

/initialize 



/s 

A 

/u 

/v 

% 



/, 



/tab 



Scope routines 



±5.1 






400000 


/• 




501004 


217200 


/J 




3762 02 


421000 


A 




1774 





A 




370047 


401170 


/m 




370040 


236000 


/n 




161044 


216000 


/o 




776IIO 
414110 


44l4oo 


$ 




477640 




370040 


202000 


h 




add 


jmp 1 cd 






add 


jmp i cd 






20100 


402010 


/- 




1012 


107000 


/) 




2010 


40201 


/; 




342 


120200 


/( 




add 


jmp i cd 






161044 


237100 


/a 




377-L04 


414000 


/b 




xbu.044 


£1L£UU0 


/c 




lb'1044 


237500 


/d 




lbl'245 


2*^6000 


/e 




'^1760 


440400 


/f 




415114 


457bOO 


/g 




376100 


434000 


/h 




7 


500000 


/I 




add 


jmp cdc 






4 





/. 




add 


jmp cdc 1 






add 


jmp cbd 


/backspace 




add 


jmp cd8 2 


/c.r. with no 


line feed 


add 


jmp cd8 







Scope routines 



16.1 



add 

30 



4010 
104422 

70402 

70020 

20242 

1012 

10027 

add 

add 

51253 

add 

add 

add 

20102 

4015 

115114 

2017 

177004 

764 04 

177003 

306240 

6047 

303214 

add 

5C241 

add 

add 

add 

add 



jmp c7d 

600 

140000 

40401 

110434 

10034 
100434 
120200 
105010 
740404 
jmp i cd 
Jmp i cd 
305010 
jmp i cd 
jmp ini 
jmp i cd 

6O70IO 
42206 
462262 
740201 
20077 
10037 
20077 
405143 
401003 
461303 
jmp i cd 
205024 
jmp i cd 
jmp i cd 
jmp cd7 
jmp 1 cd 



/space 
/" 

r 

/v 

/A 

/< 
/> 
/t 

/stop 
/pointer 

/initialize 



% 



/= 



/tab 



17.1 
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601004 


20100 


/- 


101004 


20077 


G 


376101 


210501 


/K 


377004 


20100 


/? 


376020 


200577 


M 


376o4o 


4 04177 


/N 


175014 


60276 


/o 


376110 


442206 


/? 


175015 


50336 


A 


376111 


452306 


M 


add 


jmp i cd 




add 


jmp 1 cd 




20103 


702010 


A 


1014 


77600 


/] 


1774 


740000 


X 


60200 


A 


add 


jmp i cd 




370221 


44574 


A 


377114 


462266 


/£ 


175014 


60242 


/? 


377014 


60276 


/P 


377114 


462301 


/? 


376110 


442201 


/? 


175015 


64262 


/? 


37 61 00 


402177 


/H 


1017 


76C200 


/I 


add 


jmp cdc 




lo424o 


405042 


/x 


add 


jmp cdc 1 




add 


jmp c6d 


/ba 


add 


jmp cd8 2 


/c. 


add 


jmp cd8 





/backspace 

with no line feed 



Prime tester 
5 February ±972 



1.1 

Prime tester 

/prime tester 

/nam or lam 

/preserves index register 

/call with jdp prime, and an odd positiOe number in ac 

/if number is prime, routine skips 

./if number is composite, returns with smallest divisor in location d, 
Vzero in io, (number /d)-d in ac, and does not skip 

/I is prime 

prime, 
c&c n 
law 1 
c&c id" 
a, 
law 2 

adm d 
law 1 

mul n 

div d 

hit 

sub d 

sniVsma 

jmp i prime /composite 

szm 

jmp a 

idx prime 

jmp i prime 
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/square root 

/jdp sqrt 

/input in Al, binary point to right of A(0 

/output in A, binary point to right of A(0 r 

/(alternatively, input point to right of I {16), 
/ output point to right of A(17)) 



sqrt, 



sql, 



sq2, 



srt, 





spa 
ZAI 
AVIP 



jmp 


1 sqrt 


dac 


ti 


dio 


12 


ZXP 




sxx 




rcl 


Is 


spa 




jmp 


sq2 


rcl 


Is 


sma 




jmp 


sql 


rcr 


Is 


rcr 


Is 


scr 


Is 


add 


(174000 


cxx 


! 


sar 


Is 


sxx> 


jmp 


.-2 


CXX 




sar 


Is 


dac 


t3 


lac 


tl 


lio 


t2 


scr 


2s 


div 


t3 


hit 




add 


t3 


SXX> 


jmp 


srt 


jmp 


i sqrt 



/initial guesss = x/2 + .47 



/Ixr (-1 



X 



X 



^2_ 



1 - I 
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/square root of sum of squares 



rsq, 









sar 


1 




sir 


1 




dac 


rss 




lai 






spa 






cma 






dac 


rst 




mul 


rss 




sub 


rss 




spa 






cma 






add 


rst 




dac 


rsu 




lac 


rss 




mul 


rss 




dac 


rss 




dio 


rsv 




lac 


rst 




mul 


rst 




a dm 


rss 




swp 






adm 


rsv 




lia 






and 


(i 




adm 


rss 




div 


rsu 




nop 






add 


rsu 




cli 






rcr 


1 




dac 


rsu 




lac 


rss 




lio 


rsv 




div 


rsu 




nop 






add 


rsu 




cli 






rcr 


1 




dac 


rsu 




lac 


rss 




lio 


rsv 




div 


rsu 




nop 






add 


rsu 




jmp 


i rsq 


rss, 







rst, 







rsu, 







rsv, 







constant 






start 







/first guess - about 4 bits 



/one iteration, Newton 1 s method 



/second iteration - scope accuracy 
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/sin and cos 

/J ery accurate 

/On return, A = sin, I = cos 



yOjAS 



IV\ 



\olSS\ otaxw^ 



sin, 



blah, 



xl, 
x2, 

x3, 

consta 
start 





add 

cli 

rcl 

T3X 

rcr 

sub 

dac 

mul 

dac 

scl 

dac 

sar 

sub 

mul 

mul 

add 

dac 

3ac 

sar 

sub 

mul 

mul 

scl 

add 

mul 

scl 

lio 

xct 

jmp 



(100000 

2s 

2s 

(100000 
xl 
xl 
x2 

3s 

x3 
4s 

(240573 
(373257 
x3 

(377777 
x3 
x2 

2s 

(205044 

(2370IO 

x2 

2s 
(311040 

xl 
2s 
x3 

i blah 
i sin 



nop 

craaVswp 
cmaVcmi 
cmiVswp 
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/sine routine 

/jdp sin with argument in AC 

/Very accurate 

/fast - max 227 usee 

/on input, pi=3 77777 

/on output, 1=377777 



sin, 



sc2, 



sc3, 



xl, 
x2, 

cos 








add 


(100000 


TAA3X 


and 


177777 

(100000 


sub 


dac 


xl 


X+I<M 


jmp 


sc2 


mul xl 


scl 3s 


dac 


x2 


sar 


4s 


sub 


(240573 
(373257 


mul 


mul x2 


add 


(377777 


jmp 


sc3 


mul xl 


dac 


x2 


sar 


2s 


sub 


(205044 
(237010 


mul 


mul x2 


scl 


2s 


add 


(311040 


mul xl 


scl 


2s 


TX>P 


cma 




jmp 


1 sin 

















add 


(200000 


jdp 


sin 


jmp 


i cos 



/testing X'R bit 1 
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/natural log 

/integer input, output point after bit 4 



log, 



lob, 



lot, 








spq 




hit 




ZXP 


/or ZUCP 


SXX 




rcl Is 




sraa 




jmp .-3 




rcr Is 




X-*AX 




mul (-13056 


/In (2) 


scr Is 




dio t2 




TXA 




dac tl 




ral 6s 




and (16 




TAX 




lac i lot 




a dm t2 




lac tl 




mul i lot+1 




rcl 3s 




sma 




jmp lob 




sub (400000 




sar 4s 




add t2 




add (274420 


A7.X In {2) 


jmp i log 




-ulx6£4 


730UO 


-1004^ 


65000 


-b37^ 


60000 


-5 061 


54000 


-3444 


50000 


-2600 


46000 


-1336 


43000 


-374 


41000 
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Double precision routines 

/double precis unsigned square root, +0 to +777777 777777 
/binary point to right of io bit 17 > jda sqt 
/unsigned answer in ac with binary to right of bit 17 
/if binary point shifted 2 bits in input, shift it 1 in ans 



sqt, 







dap sqx 




dio sq4 




lio sqt 




law i 22 




dac sq3 




dzm sqt 




dzm sq5 


sql, 


sad (-11 




lio sq4 




lac sq5 




ral Is 




dac sq5 




lac sqt 




rcl 2s 




sza i 




jmp sq2 




dac sqt 




sub sq5 




sub sq5 




sub (1 




spa 




jmp sq2 




dac sqt 




idx sq5 


sq2, 


isp sq3 




jmp sql 




lac sq5 


sqx, 


jmp . 


sq3, 


-22 


sq4, 





sq5, 






2.1 
Double precision routines 

constants 
start 
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/DECIMAL PRINT INTEGER. 
M. Alter, 6/9/64. 

/Enter with number in AC, jda dpi. 
/Assumes signed integer, binary point to right 
/The + sign, and all leading zeroes except one 
/ the decimal point, print as space. 

/Prints sign and 6 digits, no decimal point. 
/Exits with AC and 10 destroyed. 



of bit 17. 
just to the 



left of 



dpi, 









dap 


dix 




stf 


1 




law 


ddv 




dap 


di4 




cli 






lac 


dpi 




sma 






jmp 


dil 




cma 






lio 


(char r-J 




dac 


dpi 


dil, 


tyo 
cla 






lio 


dpi 




ril 


Is 


di4, 


div 
hit 


ddv 




dio 


dpi 




lia 






idx 


di4 




sas 


(div ddv +6) 




jmp 


di2 




elf 


1 




sni 






lio 


ddv+1 




tyo 




dlx, 


jmp 





di2, 


sni 






jmp 


di3 




elf 


1 




jmp 


dil 


di3, 


szf 


1 i 




lio 


ddv+1 




jmp 


dil 


decimal 






ddv, 


100000 


100000 


1 




octal 






constants 






start 







10000 



1000 



100 



» 
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/DECIMAL PRINT FRACTION. 

/R. Alter, 6/9/64. 

/Enter with number in AC, jda dpf, lac (n). 

/Assumes signed fraction, binary point to right of bit 0. 

/If n is positive, prints sign, decimal point, and |n| digits. 

/ If n is negative, does the same except the sign is 

/ completely suppressed. |n| may not exceed 6. 

/The + sign prints as space. 

/Maximum error = 2 in the last digit position printed. 

/Exits with AC and 10 destroyed. 



dpf , 









dap 


dfx 




xct 


i dfx 




stf 


1 




spa 






jmp 


.+3 




elf 


1 




cma 






sub 


(1) 




dac 


dfn 




add 


(7) 




srna 






jmp 


.+3 




law 


i 7 




dac 


dfn 




idx 


dfx 




cli 






lac 


dpf 




sma 






jmp 


dfl 




lio 


(char r-) 




cma 






dac 


dpf 


dfl, 


szf 
tyo 


1 i 




elf 


1 




lio 


(char r. ) 


df2, 


tyo 






isp 


dfn 




jmp 


.+2 


dfx, 


jmp 







lac 


dpf 




mul 


(12) 




rir 


Is 




dio 


dpf 




sza 


i 




law 


char rO 




lia 






jmp 


df2 


dfn, 







constants 






start 
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Double precision routines 

/DECIMAL PRINT DOUBIE. 

/Enter with number in AC and 10 in mpy format, jda dpd, lac (n), 

/ Bit 17 of 10 is ignored. 

/Assumes signed number, binary point to right of bit 17 of AC. 

/Prints sign, 6 digits, decimal point, and |n| digits. |n| may 

/ not exceed 6. Sign of n is ignored. 

/The + sign, and all leading zeroes up to but not including the 

/ one immediately to the left of the decimal point, print 

/ as space. 

/Maximum error = 2 in the last digit position printed. 

/Exits with AC and 10 destroyed. 

dpd, 

dap ddx 

lac dpd 

ral Is 

rcr Is 

lac dpd 

dio dpd 

jda dpi 

xct i ddx 

sma 

cma 

dac ddl 

lac dpd 

jda dpf 

lac ddl 

idx dds 
ddx, jmp 

ddl, 

start 
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double decimal print* rwg 
ddd, repeat 10., 69 



ddp, 











dap 


ddx 






law 


ddd-1 






dap 


dpp 






lac 


ddp 






swp 








spi 








cma 








dac 


ddl 






lio 


(charac r- 




ddl, 


idx 


dpp 






lac 


ddp 






spa 






dpo, 


tyo 
spa 

croa 








mul 


(1 






div 


(10. 




ddl, 











dac 


ddp 






lac 


ddl 






rcr 


9s 






rcr 


8s 






div 


(20. 




dd6, 


dio 

swp 


ddd-1 


/constant 




sub 


(10. 






xor 


(400000 






rcl 


1 






sar 


1 






spa 








add 


(10. 






dio 


ddl 






swp 








ior 


ddp 




dpp, 


dio 
sza 


• 






Jmp 


ddl 




ddo, 


lio 

sni 


i dpp 






lio 


(charac rO 






xct 


dpo 






law 


i 1 






add 


dpp 






dap 


dpp 






sas 


dd6 






jmp 


ddo 




ddx, 


jmp 


• 
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{k £le^<J"& 



OATING POINT PACKAGE 




Several programs are available to handle floating point data, 
of the ones mentioned here are in t h e file cabinet and on 



pfiicrotap -e- 4, under the name 'newfloat'. That file has descriptive 
subfile names, ^so you can pull out the desired functions. The 
name 'newfloat' is not meant to exclude other useful 
floating point data manipulation programs. 

'Newfloat pack* includes the basic arithmetic functions and 
some macros. It must be edited into ET before the main program that 
you write to use it. Once you have done that, your program can get 
data by using s jsp fip tyl', except that e notation is not accepted. 

The datum is stored in two consecutive words. The first word consists 
of a sign bit, a nine-bit exponent, and the most significant 8 bits 
of the mantissa. The second word is the least significant 18 bits of the 
mantissa. The exponent field of the first word contains 400 more 
than the actual exponent, A negative number Is represented by the 
36 bit one's complement of the corresponding positive number. 
To type the datum back out, use *jda fop tyo\ 

These input-output instructions 

have been made into macros. Writing ^^^tech x x will accept 
the typewriter, echo It, carriage return^TfuT^Teave the value 
and In (x,x+7). Writing 

current contents of AI, which are lost, 
ar i thmetic^operations may proceed. To put 
AI, write s load_x_l. To add the variable 



a value 
In AI 



from 



type 



the 
in, 
nto 



'jdp fad 

x s or 'mulby X s . To divide, write 



x N or 'd iadd 



To mu I tl p I y by 



will type out 
mce the data Is 
the va lue of (x,x+7 ) 
(x,x+7) to the 
contents of AI, write 
(x,x+7), write j dp fmp 
'jdp fdv^ x x or ' divby x* . To subtract, "£syja., x? , > When done, 

tore x puts the value Into (x,x+7). [Don't forget to dimension each 
varlab le with 2 spaces — nor to double array dimensions.] To 
negate the contents of AI, write 'cma^cm!', or 'jiegate/ . 
To make subscripted reference to arrays, the Index register may 
be used to index 'load' and 'store', That is, 'load I x' wi I I 
cause the address to be indexed. This will not work for the other 
macros. To index these, the effective address must be computed 
and placed under the call to the appropriate subroutine. 
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FUNCTIONS AVAILABLE 

To fix a number, put it In the AI and write N jdp fix 0* . 
Change the la tter^constant If you want a scale factor. To float 
an integer in A, *jdp f lo\ 

There are some useful constants predefined, such as x zero, 
one, two, three, four /ten, pi , 2pl s , etc. 



(exp) 
(log) 
(ran) 
(cossi n) 
(cossin) 
(sqrt) 
(besse I ) 
(besse I ) 
(besse I ) 



s jdp 

*jdp 
*jdp 
*jdp 
*jdp 
jdp fsr 



exp 
ln\ 
ran* 
cos* 
si n K 



To use the following functions, the appropriate newfloat file 
must be appended after the pack. The programs give further 
(self) description. 

will replace the value of AI by e to that power, 
will replace the value of AI by Its natural log. 
replaces by a random number In (-7.0, +7.0). 
replaces the (radian) value by Its cosine, 
replaces by the sine, 
takes the square root 
*jdp j0* replaces by Jo(AI), zeroth order Bessel fen. 
*jdp j7* Bessel of first order (and first kind) 
s jdp j2 x of second order 
(gamma) v jdp gam* replaces by the gamma (factorial) fen 
(phi) *jdp phi* replaces by std. normal distribution 
(atan) jdp^atn replaces by the arctangent. 
(shellsort) *jdp sortac' with (A) = base address, (I) = number 

^ wiLl^sort a vector of floating numbers, 

(ati) jdp ixp arg computes Altarg, arg integer. 

(atb,exp, log) s jdp fxp arg* computes Altarg, arg floating. 

W x2 y2' with x7 In the AI will draw 

from (x7,y7) to (x2,y2), using + and - 7.0 
scope limits, a nd truncating everything outside. 
? w I n 
"OS 'dj^S^to x^y'^Denuo'^^nd ~ Dje-tfdwrT telecontrol rl&rffcing 



'jdp fds 



a line 
as the 




dj^sftto x,y^ ?, penup'^a'nd 'p^eilaiwn' tp^control 
e- y arc d^rrihprkin ^owf Inaf pl^ lpark" 



There is also an ordinary differential equation integrator under 
the name newfloat ktm . The program includes a description of how 
to call It. It features automatic step-size control, using a Kutta- 
Merson scheme. 

Everything preserves your index register and flags. 



V 



x 







>V * ** irk ** icic -kit irk irk irk irk Vo'c ** *** * ** irk irk -kit : 



INDEX MODE before calling any 



The calling program must enter 
of these subroutines. 

■k-kirk-k irk irk irk **** ** ** **** irk -kic -kiricic -kic irk -k 

Feel free to add functions, or to delete parts in your own copy, 
but please DO NOT MODIFY OR DELETE THE ONES MENTIONED HERE. 

After all, you do not need to use ones you consider inappropriate. 
An example program is appended, 

2 Feb 7 96 9 

work done by William Ackerman, Rory Thompson, and Charles Landau, 
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/example program 

/examine amplitude and phase shift of a simple difference scheme 

dimension x(2 ),y( 2 ),h(2 ), d lx( 2 ), d I y(2 ) 

beg/ 




getech h 



bg7, load one 
store x 
load zer 
store y 

bg2, load x 
mulby h 
store dly 
load y 
mulby h 
negate 
store dlx 

bg3, load y 

floadd dly 
store y 
load x 
floadd dlx 
store x 

bg4, jdp fds 

y 
x 

y 

szs 10 

jmp beg 

jmp bg2 

constants 
var iab les 
start beg 
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/floating pack 3 june 7969 

/norma li ze 
/79A+20N usee, N = number of steps 

f nm, 








dac 


fmp 


X-AX 




dac 


f xr 


lac 


f nm 


dac 


fad 


TXXA<M 


cm a 




sar 


Ss 


sub 


(7 


X-A* 




spa 




cma* 


cmi 


sc I 


9s 


A+IP 


•I 


2XP 




SXX<M 


jmp 


. + 5 


re I 


7s 


sma 




jmp 


.-4 


rcr 


7s 


CXX 




scr 


9s 


X^AX 




dac 


ft7 


jmp 


fdS 



/-exp-7 



/fdv will need this 



6 J 
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/fix 






/795 


usee 


fix, 









dio 


fmp 




TAAI>P 


fx7, 


cma 






sar 


$s 




aam 






sub 


fix 




sub 


(427 




szm 






cla 






add 


(43 




spa 






c la 






rar 


3s 




add 


(add fx2+2 




dap 


fx2 + 7 




and 


fx7 




ral 


3s 




add 


(fa4+2 




dap 


fx2 




idx 


fix 




lai 






lio 


fmp 




sc I 


9s 


fx2, 


xct 


• 




jmp 


• 




repeat 4, scr 




jmp 


i fix 


/f Loi 


at 




/540- 


-20*| 


[log2(AC)J 


flo ( 


* 

lia 


- 




si r 


9s 




scr 


9s 




xor 


(210400 




jdp 


f nm 




jmp 


I flo 


f xr, 







fts. 







ft6, 







ft7, 







fts, 







ft9, 








Ss 



usee 



7 1 
Floating point package * 



steps to normalize result 



/mu I tip I: 


V 


/665 


+ 20N 


usee 


/N - 


number of 


fmp, 









dac 


ft* 




dio 


ft9 




X-AK 




cmi 






sir 


Ss 




dac 


f xr 




aam 






Ixr 


fmp 




idx 


fmp 




dac 


fad 




lac 


i 




sma 






cma 






sar 


Ss 




add 


(377 




A+IA 




dac 


ft7 




lac 


ft* 




lio 


ft9 




sc I 


9s 




scr 


5s 




ri r 


7s 




dac 


ftS 




dio 


ft9 




lac 


i 




lio 


i 1 




scl 


9s 




scr 


5s 




r i r 


7s 




dac 


ft6 




dio 


ft5 




mu I 


fU 




sc I 


7s 




sa I 


Ss 




dac 


f nm 




ri r 


2s 




dio 


fmp 




lac 


ft* 




mul 


ft9 




adm 


fmp 




r i r 


4s 




dio 


ft4 




lac 


ft5 




mu I 


ft* 




adm 


fmp 




ri r 


4s 




swp 






adm 


ft4 




lac 


ft5 




mu I 


ft9 




scr 


3s 




add 


ft4 




TAAX 


i 




scr 


6s 



/-exp(AC) 



/-exp(mem) 



/save sign of ac 
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scr Ss 

add fmp 

A$X<M 

jmp ,+4 

TAAX>P 

CI 1 1 

cma /change +0 to -0 

scr 8s 

add fnm 

A$X<M 

jmp .+4 

spa 

CI I | 

cma 

sc I 7s 

jmp fnr 



Floating point package *•* 

/add 

/665+20U usee 

/N = number of steps to normalize result 

fad, 

dio fU+7 
dac ft* 
X-AI 
dac fxr 
aam 

Ixr fad 
lac i 
sma 
cma 

sar 8s 
A-*IA<M 
cma 

sar 8s 
dac ft7 
AMIA> 
jmp . + 3 
dio ft7 
CAA| 

lxr <ftl /must be 75 bit address 

add (47 

spa 

c la 

scr 3s 

add (fa5 

dap fa3+7 

cla 

re I 3s 

add (fa4 

dap fa3 

law I 2 

adm ft7 /-exp-7 

lac I 

lio I 7 

scl 9s 

scr 3s 

and (377777 

sc I 7s 

r *i" fs /each register begins with zero 

dlo ft6 
dac fnm 

lac (ft« /IS bits 

aam 

xor fad 
A$XX 
lac I 
lio i 7 
sc I 9s 
fa3, xct . 
jmp # 

fa4, fan_=9s r epeat 72, scr fan fan_=fan>2 

fa5, repeat 4, scr 8s 
and (377777 
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fd3, TAAX /right half 



fnr, dac fmp /save sign 

/-exp-7 (must be _< -7) 

/exponent underflow 



/fraction is zero 



a dm 


f nm 


idx 


fad 


cla4 


-swp 


rcr 


H 


lio 


f nm 


add 


ft6 


TAAX 




and 


(377777 


A->IA>P 


SII 




HX> 


t 


and 


(377777 


TX<M 


jmp 


. + 3 


SAA>P 


SII 




ri I 


7s 


re I 


7s 


scr 


2s 


scl 


2s 


dac 


fmp 


Ixr 


ft7 


SX_< 




jmp 


fuv 


spa 




cma<l 


'Cml 


AUF 


1 


2XP 




SXX<M 


jmp 


. + 5 


re I 


Ts 


sma 




jmp 


.-4 


rcr 


/s 


CXX 




scr 


Ss 


sc I 


7s 


r i r 


1s 


SII>P 


SAA 




ri I 


Is 


scr 


2s 


dac 


fdv 


ral 


9s 


X-*A<M 


jmp 


• + 5 


SAX 




lac 


fdv 


scr 


7s 


X-AX 


L 


Ixr 


ral 


Ss 


A+XA>P 


jmp 


fov 


Ixr 


fmp 


TX>P 


> 


cma* 


cmi 


Ixr 


f xr 


jmp 


i fad 



/check sign of left half 

/carry into right 

/in case that turned on sign of right half 



/carry into left half 
/wrap around to right half 



/put exp in XR (_> -0 ) 



/fraction overflow 
/shift fraction down 



fd£, r a .L.^? /exponent 

/overf low 
Ixr fmp 
TX>P 

fnx, 
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fov, lac (377777 
c I ?*cmi 
jmp fnx-3 

fuv, cla+cll 

jmp fnx-3 
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/d i vl de 
/S35+20N usee 
/N = number of 



steps to normalize result 



fdv, 



ft4. 





dac 

d io 

TXA 

aam 

Ixr 

I io 

Ixr 

X-AX 

jdp fnm 



ft5 
ft6 



fdv 
i 7 
? 



sc I 
dac 
idx 
dac 
r I r 
dio 
lac 
sma 
cma 
sar 
sub 
adm 
lac 
lio 
sc I 
scr 
dfv 
jmp 
dac 
swp 
mu I 
diy 


dac 
lac 
mu I 
div 
jmp 
mul 
add 
lio 
sir 
Ixr 
I«-XI 
spa 
I+XI 
dio ft5 
lio fnm 
si I Ss 
si I Ss 
A*-XA 
I«-XI 
A + IA 
lio ft5 
jmp fd3 



9s 

ft9 

fdv 

fad 

7s 

ft* 

f t5 



Ss 

(402 

ft7 

ft5 

ft6 

9s 

7s 

ft9 

fov 

fnm 

fc7 
ft9 

ft5 

ft* 

(-200000 

ft9 

fov 

fnm 

ft5 

fnm 

7s 

(377777 



/division by zero 



1200000 
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/floating input 



f ip, 



fi^ 



f!2, 



dap 
TXA 
dac 
law 
TAX 
dap 
d2m 
dim 
law 
dac 
xct 
lac 
A^I 

jmp 
sad 
cla 
sad 
jmp 
sas 
jmp 
law 
dap 
jmp 
add 
dac 
and 
sas 
jmp 
lac 

jdp 

ten 

jdp 

f 17 

dac 

dio 

law 



jmp 



fi7 

f lo 
i 7 

f im 
f \8 
f\9 
i-£ 
f ?2 

fiS 
AP | 
f \4 

(to 

(73 

f\5 

(54 

. + 4 

600 

fim 

fil 

(204000 

f'\7 

(-360 

f\7 

f\4 

f 19 

fmp 

fad 

f I 9 

f\S 
i 7 

fi7 



10 
/. 

/- 



v- A 



fi4, idx fi7 
lac fi9 
SXX<M 
jmp fim 
jdp fdv 
ten 
jmp .-4 

fim, skp i 

cma^cmi 
Ixr flo 
jmp i fi7 
fi5, law i-A+XX 
jmp fi7-7 



//*? - \aA 
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/floating output 

/number of digits printed can be changed at fop+6 and foS+4 



fop, 









dap 


fo9 




SAA 






dap 


fox 




law 


i-2 




dac 


fig 




law 


i 8. 




dac 


f I 7 




lac 


fop 




dzm 


f lo 




A<-'IM| 




2AI 






d io 


f \8 




sma 






jmp 


fo7 




cma^cmi 




dio 


f\8 




lio 


(charac r- 




xct 


fo9 




lio 


f IS 


fo7, 


dac 


f 19 




jdp 


fix 











lio 


(jdp fmp 




sza 


I 




jmp 


.+4 




law 


i-CAA 




lio 


(jdp fdv 




dac 


fig 




dio 


fo2-7 




law 


fo5-2 




dac 


fo2 




law 


100 


fo2-< 


$, 






law 


2 




a dm 


fo2 




lac 


f i9 




lio 


f \8 









fo2, 









dac 


ft9 




jdp 


fix 











xct 


fig 




TA<P 




jmp 


f33 




lac 


f t9 




lio 


fmp 




dac 


f i9 




dio 


f \8 




AIIP 


fo3, 


law 




fig, 









a dm 


f lo 


f33, 


xct 


fo3 




rar 


7s 



/number of digi ts + / 



/minus zero 



/number is < 7 
/number Is > 7 



dap fo3 



/jdp fmp or fdv 
/some power of ten 

/this contains a dio fmp 



/don't change exponent If exactly zero 
/2, or CAA if num was _> 7 
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jmp 


fo2-6 










xct 


fig 










spa 












jmp 


fOO 










Law 


i i 




/ni 


jm 




adm 


f lo 










idx 


f 17 










lac 


f f 9 








I 


io m 










jmp 


fo7 + 6 








fOO 


, llo 
xct 

jmp 


(charac i 

fo9 

foS 


m 






/f 1 


9, fli 


I are now 


stri 


ct ly 


< 



/number was _> / 



fo7, 



fo9, 



foS, 



fox, 



lac 


f 19 


llo 


f 1* 


jdp 


fmp 


ten 




dac 


f 19 


dlo 


fl* 


j dp 


fix 







TAI 


1- 


lio 


(charac 


xct 


• 


rl r 


5s 


law 


202 


re l 


9s 


cma<l 


'C U Icml 


jdp 


fad 


f i9 




dac 


f 19 


law 


7 


AilA 


dac 


f IS 


i dx 


f lo 


sza 


i 


jmp 


fOO 


i sp 


f i7 


jmp 


fo7 


lac 


f lo 


CAAJ 


M< 


add 


(7 


TAA= 




A$I>P 


jmp 


fox 


cii 




xct 


fo9 


lio 


(charac 


xct 


fo9 


cli 




xct 


fo9 


lio 


(charac 


CAA>P 


CAA| 




xct 


fo9 


dac 


ft9 



rO 



/compensate for truncation 



/number of di gi ts 



re 



r- 
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dim 


ft* 


dpp, 


dac 


ft7 




mu L 


(7 




di v 


. + 7 




12 






sas 


fts 




jmp 


dpp 




sni 






lio 


(charac 




xct 


fo9 




lac 


ft7 




dac 


ftS 




lac 


f t9 




sas 


ftS 




jmp 


dpp 




jmp 


fox 


fo5, 


352702 




265635 




233276 




275676 




207234 


nun, 


203710 


ten/ 


202240 


f 19/ 







f I J, 







f*7. 








r0 



360175 110*64 

613267 110*32 

067447 I10M6 

70200 /10tS 

fc7 ,200000 
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/macros for new pack, 27 aprll 68 

define getech a 

Jsp flp 
tyl 
dac a 
dlo a+7 
jda fop 
tyo 

c 1 1 1 cm I 
tyo 
Lac a 
llo a+1 
terml nate 

define load a 

lac a 
llo a+J 
terml nate 

define fadd a 

j dp fad 
a 
terminate 

equals floadd,fadd 

define mulby a 

jdp fmp 
a 
terml nate 

define store a 

dac a 
dio a+7 
terminate 

define fsub a 

negate 
jdp fad 
a 

negate 
terml nate 



def i ne 

j da fop 

tyo 

terminate 



type 



define di vby a 

jdp fdv 
a 
terminate 

negate=cma+cmi 
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def ine 

c 1 1 Icmi 

tyo 

terminate 

def ine 

jsp txx 
text a 
lio (236 
tyo 
terminate 



carrig 



a 



\>5o /df 



write a 



txx, dap txy 
a am 

lio txy 
idx txy 
lac (607600 
rcl 6s 
sad (lai 

txy, jmp 

sad . 1 

jmp txx / 

swp 

tyo 

lia 

jmp txy-3 
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/constants 






zero, 








2er=2ero 






one, 200600 







two, 201200 







three, 


201300 





thr = three 






four, 


201600 





for=four 






five, 


201640 





flv=five 






/ten and hun (100.) 


exist elsewh 


ere 


pi, 201311 


37553 




pi2, 200711 


37553 




2pi, 201711 


37553 




m/, 577177 


777777 




m2, 576577 


7 77777 




haf, 200200 







pt7, 7767/4 


637463 


10.1 


P07, 175243 


656057 


I0 % 01 


rt2, 200665 


77775 


/sq root of 



2.0 

dimension tem(2 ), tmp(2 ) /so always available 
dimension f b7 (2 ), fb3(2 ), f b5 (2 ) /used by functions 
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/sqrt, 4 feb 7 96 9 
/Newton-Raphson method, 
/ca I I by jdp f sr 
/sqrt(AI) -* AI 
/5.7 msec 



3 iterations 



fsr, 



f sa, 








spa 




negate 


AHP| 


jmp 


i f sr 


dac 


f \7 


dio 


f lo 


scr 


7s 


rar 


Ss 


spa 




add 


(202000 


ral 


Ss 


add 


(100172 


dac 


f 19 


law 


i 3 


dac 


fix 


dio 


f i9 + 7 


lac 


f i7 


lio 


f lo 


divt 


>y f ! 9 


floadd f!9 


sub 


(400 


dac 


f 19 


isp 


fix 


jmp 


f sa 


lac 


f 19 


jmp 


i fsr 



/approximate sqrt(x) by 
/ x/2+ t 4750 in £.5,7) 
/or X/2+.4645 in [7,2) 
/relative error < .0355 
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/log, 18 dec 7 96* 
/ca I I by jdp fin 
/ln(AI) - AI 
110 msec 

ln2, 200261 3 44137 /. 6931471 (ln(2)) 

fin, 

jdp lg2 

jdp fmp 

ln2 

j mp i fin 

/log base 2 

/Hastings, Approxs. for dig. comp., pi 66 
/cal I by jdp lg2 
/log2(AI) - AI 
110 msec 

lg2, 

jdp fnm 
spq 

hit /zero or negative 

and (377 

dio fb3 J 

Ho (-401 

sub (265 

sma 

SI 1 1 />7/sqrt 2 

add (400 /<7/sqrt 2 

add (200265 

dac fb3 

lai 

add ft7 

jdp flo 

store f b7 

lac one 

cli 

floadd fb3 

store fb5 

c li *cmi 

lac m7 

floadd fb3 

divby fb5 

store fb5 

mu I b y f b 5 

store fb3 

mulby flh 

floadd flh+2 

mulby fb3 

floadd flh+4 
mulby fb3 
floadd flh+6 
mulby fb5 
floadd fb7 
jmp i lg2 
flh, 177736 256453 / .4342597 

200223 460041 /. 5765385 

200366 161113 / .961 8007 

201270 5 24353 12.885390 
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/exp, 11 nov 7 96 9 
/ca LI by j dp exp 
/et(AI) - AI 
II msec 



exp, 



j d p f mp 
. + 3 

jdp f2x 
jmp i exp 



200670 



524355 II .4426950409 (log2(e)> 
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/2tx 

/jdp f2x 

/2t(AI) -» AI 

/Hart et a I, Computer Approx's, number 1063 

/approximation is accurate to . 7SE-9 relative 

/6.5 msec 

f2x, 

store fb3 

div (204400 

hit /| exp| _> 256. 

load fb3 

jdp fix 

-*• 

dac fb3 

and (777400 

spa 

sub (7 
/separation of integer and fraction parts is bad by one bit 
/here if arg is negative, maybe someone would like to fix it 

dac f b 7 

law 377 

and fb3 

xor fc7 1100000 

j dp f nm 

floadd f2a+72 
/if fraction part = 0.5, will get zero here, divide overflow 
/will occur, seems to get correct result exactly 

store fb3 

load f2a 

divby fb3 

floadd fb3 

store fb5 

load f2a+2 

divby fb5 

floadd fb3 

floadd f2a+4 

store fb5 

load f2a+6 

divby fb5 

floadd f2a+70 
add fb7 
jmp i f2x 

f2a, 202646 432240 / 20.81 $923793 

203720 400752 / 104.25093406 

575165 2 47724 1-17 .334004949 

574473 706557 I -49 .027969777 

577112 7 66062 /-I .41 421 35623 (-sq rt of 2) 

577577 777777 /-0.5 
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/cossin, 16 feb 7 96 9 

/jdp sin or jdp cos 

/sin(AI) or cos(AI) ■* AI 

/Hart et al, Computer Approx's 

/approx no. 31 SO (7th degree poly) for sin 

/approx no. 3700 (6th degree poly) for cos 

/both are accurate to .9E-* (relative) 

II \9 msec 

cos, 

floadd pi 2 
dac fb7 
lac cos 
dac sin 
lac fb7 
jmp sin+7 

sin, 

di vby p i 2 

store fb3 

jdp fix 



TXI 

dio fb7 

CAAKM 

SAA| 

sub (7 

sar 7s 

scr 7 s 

A$IX 

sc I 2s 

jdp flo 

floadd fb3 
/in first or fourth quadrant 
/-I _< AI _< 1, want sin (AI*pi/2) 
/XR has sign 

A$XX 

spa 

negate 

dac fb5 /now in [0,7] 

sub (200222 /.56S8 

spa 

jmp sn2 
/use cos in [0, .4372*pl /23 

lac fb5 

floadd m7 

dac fb5 

law slt+70 

jmp sn2+7 

/use sin in [0, .56 SS*p i/2] 
sn2, law sit 

dac sic 

add (2 

dac sia 

lac fb5 

dio fb5+7 

mulby fb5 
store fb3 
jdp fmp 
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sic, 

jmp si a-/ 
sib, lac cos 

mulby fb3 

jdp fad 
sia, 

dac cos 

law 2 

a dm sia 

sad (sit+20 

jmp scg /end of cos 

sas (s!t+/0 

jmp sib 

lac cos /end of sin 

mulby fb5 

jmp ,+2 
scg, lac cos 

TX>P 

negate 

Ixr fbl 

jmp i sin 

sit, 603151 767434 /-. 00457 S1 3997 

[76643 125500 / .07967150287 

577532 504277 /-. 6459627091 6 

^00711 37552 / 1.5707963124 

602127 721431 /-. 02051 88 B6 93 

177601 667174 / .25362870276 

577142 054232 /-7 .2336989859 

200377 777777 / .99999999043 
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/elliptic functions, 1$ feb 7969 

/argument is parameter, or modulust2 

/Hart et al, Computer Approx's, numbers 7303 and 7403 

/both are form P( 7 -x )- log( 7-x )Q( 7 -x ), P and Q quintic 

/accurate to .35E-9 absolute 

/jdp elpe or jdp elpk 

/jdp elpec or jdp elpkc for complementary functions 

IIQ msec 

elpk, o 

dac el7 
lac elpk 
dac e Ipe 
law ept+24. 
jmp eiO 

e I pe , q 

dac e 17 

law ept 
el0, dac ep5 

lac e 17 

negate 

floadd one 
e 12, store el7 

jdp lg2 

store fb3 
ep3, law 2 

add ep5 

dac ep6 

load el7 

jdp fmp 
ep5, 

jmp ep6-7 
ep4, lac elpk 

mu Iby e 1 7 

jdp fad 
ep6, 

dac elpk 
law 2 
a dm ep6 
sas (ept+72, 
sad (ept+36. 

jm 5 fP ; /done with first polynomial 

sad (ept+24. 

jmp ,+3 

sas (ept+4S # 

jmp ep4 

lac elpk /done with second polynomial 

floadd fb3 

jmp i elpe 

ep7, dac ep5 
lac elpk 
mulby fb3 
negate 
store fb3 
jmp ep3 
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e Ipkc, 





dac fb7 




lac elpkc 




dac e Ipe 




law ept+24. 




jmp ecO 




elpec, 





dac fb7 




lac elpec 




dac e Ipe 




law ept 




ecO, dac ep5 




lac fb7 




jmp e 12 




ept, 173701 


025376 


175367 


766326 


176232 


032777 


176604 


644471 


177261 


343714 







174772 


613277 


175773 


1 37152 


176202 


057746 


176353 


6 26576 


177742 


711674 


200600 





173647 


646422 


175324 


750217 


175777 


334213 


176307 


122541 


176661 


343532 


177661 


344140 


174731 


445301 


175724 


540176 


175702 


515632 


176201 


210732 


176705 


626014 


200661 


3 44140 



/K(7-x) 



/E(7-x> 



/ .001472793465 
I .015135576508 
I . 037610529537 
1 .064854252062 
1 . 173286242518 
/0. 

I .007652960603 
1 . 030662347457 
1 . 031761145525 
/ .057566998484 
1 . 443152874726 
/I. 

1 . 001280405152 
1 .012997660452 
1 . 031180446501 
1 . 048623413677 
1 .086642909577 
/. 346573590280 

1 . 006639801115 
1 .025962888453 
1 . 023761224858 
I .031559431628 
1 .096578619623 
II .38629436112 



el7 / 
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/arcsincos, 6 feb 7969 

/Hart et al. Computer Approx's, number 4693 

/call by jdp ars or jdp arc 

/arcsfn(AI) or arccos(AI) •* AI 

/approximation is accurate to 2.5E-9 

/about 7 msec in 1-0.5,0.51, 72 msec elsewhere 

ars, 

dac fb7 
TXA 

dac xrs 
lac fbJ 
sma 

sxx| 

negate 
sub haf 
spa 

jmp as6 

add (200200-400 
negate 
floadd haf 
TAA_> 
hit 

jdp fsr 
CXX| 
as6, add haf 

store fb5 
mulby fb5 
store f b 7 
floadd as7 
store fb3 
load as7 + 2 
divby fb3 
floadd as7+4 
floadd fb7 
store fb3 
load as7+6 
divby fb3 
floadd as7+70 
mulby fb5 
TXX<M 
jmp .+5 
add (400 
negate 
floadd pi 2 
TX| = 

negate 
Ixr xrs 
jmp i ars 



/in 10.5,1.01 
I -x/2 
/ (7-x)/2 

/arg was not in 1-1, +7] 



/now in 10,0.51 

/uses x*P( x t2)/Q(xt2) 
/P and Q are quadratic 



/if x was _> .5 , 

/use pi/2 - 2 arcsin(sqrt((7-x)/2)) 



/restore sign 
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arc, 

jdp ars 
negate 
floadd p i 2 
jmp i arc 

as/, 577075 
577575 
576176 
577021 
177775 



766502 /-7. 5157679526 

543412 /-. 50900634073 

750207 /-4.03269S6467 

242072 /-1 ,864713X142 

576334 I .49559947479 



xrs, 
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/atan, 27 mar 7 96 9 

/ca I I by jdp atn 

/arctanCAI ) ■* AI 

/Hart et al, Computer Approx's, number 5057 

/approximation is accurate to 7.05E-9 relative 

/approx 1 1 msec 

atn, 

dac fb7 

TA>P 

law cma^cmi -opr 

dap svs /save sign 

law at4 + 2 

dac at6 

lac fb7 

xct svs 

store f b 7 

dzm atB 

sub (177611 

spa 

'<** at * /| x| _< .26795 = tan(pi/72) 

sub (76 7 

spa 

•dx atB /I xl _< 7 = tan(p?/4) 

sub (557 

spa 

•dx at* /| x | _< 3.732/ = tan(5pi/72) 

law 3 

sub at£ 

jdp flo 

mu Iby pi 6 

store fb7 /amount to be added back later 

law 2 
sub atB 
SAP | 

jmp. at3 /no correction 

ra I 7s 
add (at7 
dac at5 
dac at8 
load fb7 
jdp fad 
aU, 

store f b 7 
TXA 

Ixr at5 

lio i at9-at7+7 
Ixr i at9-at7 
X-AX 

divby fb7 

negate 

jdp fad 
at5, 
at2, store f b 7 /now | AI | _< tan(pi/72) 

mulby f b 7 

store fb3 

mulby at4 

jdp fad 
at6, 
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at7, 



svs, 



at3, 



at7, 



at9, 



at4, 



pi6, 
fb7, 



store fb5 
law 2 
adm at6 
sad (at4+70 
jmp at7 
load fb3 
dlvby fb5 
jmp at6-7 
load fb7 
divby fb5 
floadd fb7 
opr 

jmp I atn 
load fbl 
jmp at2 + 2 

200735 
200223 


207600 
200652 
200600 

177636 
200640 
201300 
200600 

200206 




547535 /sqrt(3) = 1 .7 320S0 S07S7 

32357 /sqrt(3)/3 - .5773502697 90 
/0 

/4 
525253 1413 
II 

714464 /.3103S080523 
037746 /7. 2504*75062 

000111 13.0000043545 

000000 II .0000000011 

025707 /pi/6 = .5.2359*77559* 
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/ati, 20 dec 196S 

/call by jdp i xp followed by exponent (integer) 

/Altexponent -* AI 

/$ 7.7 log2(exponent ) msec 

ixp, 

store fb3 

I i o one 

dio fb7 

dzm fb7+7 

TXA 

dac fb5 

Ixr ixp 

idx Ixp 

Ixr i 

lac i 

i io (jdp fdv 

spa 

CAA| 

Ho ixb+2 /jdp fmp 

dio ixo 
ixc, scr Ts 

TAX 

sp? i 

jmp i xb 

load fb7 

' xo / /jdp fmp or fdv 

fb3 

store f b 7 
ixb, load fb3 
ixb+2^ mulby fb3 

store fb3 

TXAP 

jmp ixc 

load f b7 

Ixr fb5 

jmp i ixp 
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/atb, 10 dec 1 96 8 

/call by jdp fxp followed by exponent 

/Altexponent ■* AI 

/base must be _> 0, 0t0 = 

/17>5 msec 

fxp, 

dac fmp 

a am 

lac fxp 

dac ,+7 

fdx fxp 

lac fmp 

S2a I 

jmp ? fxp /base is zero 

jdp lg2 

jdp fmp 



jdp f2x 

jmp I fxp 
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/cube root, 27 mar 1969 

/Hart et a I, Computer Approx's, number 560 

/followed by 3 Newton iterations 

/ca 11 by jdp cbr 

/cbrt(AI) - AI 

/72 msec 



cbr, 







store fb5 




X-AX 




dac cbx 




TXXA>P 




negate 




store fb3 




sar Ss 




SAA 




mul (7 




div ,+r 




-3 




add (7 25 




sa I Ss 




dac fbT 




law 376 




A+II 




si I Ss 




law 377 




and fb3 




AHA 




lio fb3+7 




mulby cba 




floadd cba+2 




sub fbT 




TX>P 




negate 




Ixr (-3 


ebb, 


store fb3 




mulby fb3 




store f b 7 




load fb5 




divby fb7 




negate 




floadd fb3 




divby mth 




floadd fb3 




SXXP 




jmp ebb 




Ixr cbx 




jmp i cbr 


cba, 


200232 




177750 


cbx. 





mth. 


-201300 



534046 / .60421 SI 

12057 /.4531635 

-0 



Floating point package 



l% % 1 



/ran, 16 mar 1969 
/ca 1 1 by jdp ran 
/random number in 



(0,11 



AI 



ran 








lac ;r7 
xor (311071 
add (311071 




rar 9s 




dac r7 




lac ;r2 
xor (355671 
add (355671 




rar 7s 

dac r2 




xor r7 




c I i iswplcmaicml 
re I 7s 




and m7 
floadd two 


K-100600 


jmp i ran 





29 7 
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/phi 

II apri I 68 

/ca 11 by phi 

phi =jdp ph2 



ph2. 







dac fb7 




TXA 




Jac fb3 




lac fb7 




sma 




CXXM 




cmatcml 




dlvby rt2 




store f b 7 




load a6 




mulby f b 7 




floadd a5 




mulby f b 7 




floadd a4 




mulby f b 7 




floadd a3 




mulby f b 7 




floadd a2 




mulby f b 7 




floadd a7 




mulby f b 7 




floadd one 




store fb7 




mulby f b 7 




store f b 7 




mulby f b 7 




store fb7 




mulby f b 7 




store f b 7 




mulby f b 7 




store fb7 




load haf 




divby fb7 




TXP 




jmp ph7 




negate 




floadd one 


ph7. 


Ixr fb3 




jmp i ph2 


a7, 


1766Z0 


a2, 


176255 


a3, 


775227 


a4, 


772237 


a5, 


772627 


a6, 


777264 



/if positive argument 



/fourth power 



334631 

137631 

706637 

314031 

3 36 

476631 



Floating point package 

/pull In tape to first carriage return 
/i.e., expects a title to Ignore 



30.1 



hdr, .-. 

law char rr 

arq 
nrd, hit 

rpa 

lai 

and U7 

sas (77 

jmp nrd +7 

jmp I hdr 



/use jdp hdr 

/reader not aval lab le 

/If not cr 



Floating point package 

/read one value and return it in AI 



37.7 



/use so: 

If i rst / jdp hdr 

/then, j sp fip 

hnd, .-. 

dac tmp 
rpa 
lai 

and (77) 

dac tmp+7 

sad (20 

jmp hn4 

sas (0 

j mp . + 7 
lai 

and (777 

sas (200 

jmp hn3-7 

lac tmp 

jmp i hnd 

sub (12 
spa 

j mp hn 4 

lac tmp+7 

sad (73) 

jmp hn4 

sad (54) 

jmp hn4 

sad (7 7) 

jmp hn 4 

sad (36) 

jmp hn4 

sad (7 3) 

jmp hn3-7 

jmp hnd+2 



to get in the beginning 
jdp hnd 



/so looks like tyi 

/so exit via hn4 when accept 

/0, so ok 



/exit if blank tape 

/assuming space is acceptable to fip 

/an integer, so okay 

/period 

/mi nus 

/carriage return 

/tab 

/stop code 

/ignore any other symbols 



hn4, load tmp 
jmp i hnd 

TIIX 
hn3, hit 



/change to a jump if want an exit 



Floating point package 
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/j0, Besse I function of zeroth order and first kind 
dimension rj2(2) 

dimension xv3(2 ),xsav(2 ) 
3vx=xv3 



JO, 



spa 

negate 
store xsav 
divby thr 
store xv3 
fsub one 
sma 
jmp ap94 3j 

ap947, 

mulby xv3 
store xv3 
mulby p00027 
fadd m00 394 
mulby xv 3 
fadd 904444 
mulby xv 3 
fadd m3U3S 
mulby xv3 
fadd 7p2656 
mulby xv 3 
fadd 2m2499 
mulby xv3 
fadd one 
jmp i j0 

dimension jamp(2 ), j th( 2 ) 

ap943j, 

load xsav 
jdp fsr 
store tern 
load j th 
jdp cos 
mulby j amp 
divby tern 
jmp i j0 

P00027, 
m00394, 
904444, 
m3l63S, 
1^2656, 
2m2499, 



/jdp j0 with argument in A I 

/even function 
/for ap943 



load xv3 /x .le. 3 in this approx 



/value returns in AI 



jdp ap943 /x ,g. 3 for this approx 



/phase 

/amp 

/sq rt of x 



172334 
603176 
176266 
600136 
200641 
576560 



146720 

577746 

36000 

5114 

777565 

1 



Floating point package 



33.1 



/}0 cont. 

ap943, 

load thr 
divby xsav 
store 3vx 
mu Iby 9OOOI 4 
fadd m00072 
mulby 3vx 
fadd 900137 
mulby 3vx 
fadd m00009 
mulby 3vx 
fadd m00552 
mulby 3vx 
fadd m00 77 
mulby 3vx 
fadd Q797SS 
store jamp 

a P 9432, 

mulby p000l3 
fadd m000Z9 
mulby 3vx 
fadd m00054 
mulby 3vx 
fadd 900262 
mulby 3vx 
fadd m00003 
mulby 3vx 
fadd n\04l66 
mulby 3vx 
fadd n\7S5 39 
fadd xsav 
store jth 
jmp i ap943 



/amp and phase for j and yO, 



load 3vx 



pOQOU, 

m00072, 

900137, 

m00009, 

m00S52, 

m0077, 

979711, 

900013, 

m00029, 

m00054, 

900262, 

m00003, 

m04166, 

m7SS39, 



172227 
604501 
173663 
606070 
603112 
611461 
200314 
172216 
605146 
604562 
174254 
606532 
601525 
577466 



625335 

112620 

702141 

411533 

701473 

234004 

204246 

124751 

153650 

72516 

50701 

120376 

260246 
740224 
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34.7 



/J7 

/j7, Bessel function of first order and first kind 

/requires cos and fsr 

/dimension xv3(2 ),xsav(2 ) if not using j 



store xsav 

spa 

negate 

divby thr 

store xv 3 

fsub one 

sma 

jmp ap94 6 

ap944, 

mulby xv3 
store xv 3 
mulby pOOOOl 
fadd n\00031 
mulby xv3 
fadd p00443 
mulby xv 3 
fadd m03954 
mulby xv3 
fadd pZ1093 
mulby xv3 
fadd m56Z49 
mulby xv3 
fadd haf 
mulby xsav 
jmp i j7 

pQOOOl, 
m00031, 
P00443, 
m03954, 
p27093, 
m56249^ 



/jdp jJ with argument in A I 



/out of range for ap944 



load xv 3 



/answer returns in AI 



170272 
605131 
174621 
601536 
177327 
577560 



36276 

366156 

210452 

20431 

777041 

6 



Floating point package 35,7 

/j 1 , con t 

/this can also be used for Y7, simply by using sin in place of cos below 

load thr 



ap946, 

dlvby xsav 
spa 

negate 
store 3vx 
jdp ap946b 
load xsav 
spa 

negate 
jdp fsr 
store tem 
load j th 

ap94<Sy, 

mulby jamp 
divby tem 
dac tem 
lac xsav 
sma 

jmp ap946p 
lac tem 
negate 
jmp i j 1 

ap946p, 

jmp i jT 



/approx for 3 < x 



/sqrt 
jdp cos 



/change to jdp sin for y/ 



lac tem 



Floating point package 56 • * 

/j7, cont. 



ap946b. 


• ~ • 




load 3vx 






mulby m00020 






fadd p00 7 73 






mulby 3vx 






fadd m00249 






mulby 3vx 






fadd p00077 






mulby 3vx 






fadd p07659 






mulby 3vx 






fadd p0075 






mulby 3vx 






fadd p797BS 






store j amp 






load mooo2 9 






mulby 3vx 






fadd p00079 






mulby 3vx 






fadd p00074 






mulby 3vx 






fadd m00637 






mulby 3vx 






fadd p00005 






mulby 3vx 






fadd p72499 






mulby 3vx 






fadd 2m3567 






store jth 






load xsav 






spa 






negate 






fadd j th 






store jth 






jmp i ap946b 






m00020, 


605455 


740577 


P00 7 7 3, 


773624 


757777 


m00249, 


603534 


366007 


P00077, 


772263 


267672 


P07659, 


775607 


753364 


P0075, 


766727 


302307 


mooo2 9 , 


605747 


54074 


P00079, 


773327 


207757 


P00074, 


773302 


774742 


m00637, 


603056 


765526 


P00005, 


777354 


764646 


p72499, 


776777 


773733 


2m3567, 


576557 


750750 



Floating point package 

in 



37.7 



dimension rj2(2) 
/uses recurrence: 

J2, .-. 

store rj2 
spa 
cma 

and (377 
sas (0 
jmp nzrj2 
load zer 
jmp i j2 

nzrj 2, 

J d p j 7 
divby rj2 
mulby two 
store tmp 

load rj2 
jdp j0 
negate 
fadd tmp 
jmp i j2 



j2 = 2 j 1 over r - j0 

/use jdp j2 with arg In AI 



load rj2 



Floating point package 
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/Gaussian integration subroutines 

/call either the seven pt or the nine pt 

/by jdp, with the lower limit of integration 

/in AI, and the address of the upper limit 

/in the address following the jdp. The 

/function to be integrated must be reachable 

/by xct I [second address after the jdp] 

/with the argument in AI, 

/e.g., load a jdp gauss7 b jdp fen 

/two routines are provided for Independent estimates 

dimension ml dpt(2 ), hafwd(2 ) 

dimension fcnadd(2) 



/absci ssa 
xi7, 177717 
200275 
200362 

/wei ghts 
177725 
wi7, 777703 
177617 
177204 

xi9, 777646 
200235 
200326 
200367 

177651 
wi9, 777637 
177605 
177270 
176646 



625675 
651763 
761356 



775373 
376351 
153055 
457324 

7052 
73647 
74270 
662541 

52213 

727722 

335417 

767523 

346305 



gauss7, 

store hafwd 

rpf 

nam 

dio fenadd 7 

lac gauss7 

dac midpt 

lac i midpt 

dac midpt 

SAA 

dac midpt+7 

idx gauss7 
dac fenadd 
idx gauss7 



•". /use jdp 

/lower limit of Integration 



/load i midpt to get upper limit now 



load i midpt 
fadd hafwd 
divby two 
store midpt 
fsub hafwd 
store hafwd 



Floating point package 

load zer 
store sum 



39.7 



load mfdpt 
xct i fcnadd 
mulby wl 7-2 
store sum 



/get value at midpoint 



load xi7 
negate 
mulby hafwd 
fadd mldpt 
xct i fcnadd 
mu Iby wi 7 
fadd sum 
store sum 

load xI7 
mulby hafwd 
fadd midpt 
xct i fcnadd 
mu Iby wl 7 
fadd sum 
store sum 

load xl7 + 2 
negate 
mulby hafwd 
fadd mldpt 
xct i fcnadd 
mu Iby wi 7+2 
fadd sum 
store sum 

load xi7+2 
mulby hafwd 
fadd midpt 
xct i fcnadd 
mu Iby wi7 + 2 
fadd sum 
store sum 

load xi 7 + 4 
negate 
mulby hafwd 
fadd midpt 
xct i fcnadd 
mulby wI7+4 
fadd sum 
store sum 



load xi7+4 
mulby hafwd 
fadd midpt 
xct i fcnadd 
mu Iby wl 7 + 4 
fadd sum 

mulby hafwd 
dio hafwd 



/for given length 



Floating point package 



39.2 



lio fcnadd+7 

Ipf 

lio hafwd 

jmp i gauss7 



/result in AI 



Floating point package 40 * 1 

/Gauss9 

/ca I I just as gauss7 

gauss9, .-. /use jdp 

store hafwd 
rpf 

dio fcnadd+7 

lac gauss9 

dac mldpt 

lac i midpt 

dac midpt 
SAA 

dac midpt+7 

idx gauss9 
dac fcnadd 
idx gauss9 

load i midpt 
fadd hafwd 
divby two 
store midpt 
fsub hafwd 
store hafwd 





loaa 


I zer 




store sum 




load 


I midpt 




xct 


i fcnadd 




mu lb 


• y wi 9-2 




store sum 




d2m 


;k 


gs9l, 


t 






add 


k 




add 


k 




dap 


gs9a 




dap 


gs9ap 




SAA 






dap 


gs9a+7 




dap 


gs9ap+7 




law 


w!9 




add 


k 




add 


k 




dap 


gs9b+7 




dap 


gs9bp+7 


gs9a, 


f 






negate 




mulb 


y hafwd 




fadd 


i midpt 




xct 


i fcnadd 


gs9b, 


9 






fadd 


; sum 




stor 


e sum 



law xi9 



load 



mu Iby w 19 



gs9ap, load 



Floating point package 4 * •* 

/gamma function, 24 jan 7 96 8 

/ca I L by jdp gam 

/G(AI) - AI 

/C. Hastings, Approx's for Dig, ComD. p 158 

/accurate to 2.5E7 In range [7,2] 

/approx 73,S+7.5[| z-7.5| +0.51 ms 

gam, 

store fb3 

load one 

store fb5 

load fb3 

floadd m7 
gm7, TAA>= 

jmp gm5 /too smal I 

store fb3 

floadd m7 

TAA> = 

jmp gm7 /0__<fb3<7 

store fb7 /too big 

load fb3 

mulby fb5 
gm2, store fb5 

load fb7 

jmp gm7 
gm5, load fb3 

store fb7 

floadd one 

store fb3 

load fb5 

divby f b 7 

jmp gm2 
gm7, law gaO+2 

dac gmS+4 

load ga0 

jmp .+2 
gmS, lac fad 

mulby fb3 

jdp fad 
gmS+4, Q 

dac fad 

law 2 

a dm gm£ + 4 

sas (ga0+22 

jmp gmS 

lac fad 

mulby fb5 

jmp i gam 

ga °' llttll 7 25274 I .035868343 

tO, * 71 * 47537 /-,793527m 

177766 705533 / ,482199394 

577476 227220 /-. 756704078 

100353 037153 I .918206857 

577432 264747 /-. 897056937 

200374 766077 / ,988205891 

577554 772253 /-, 577797652 

200600 11.0 



Floating point package 

mu Iby hafwd 
fadd midpt 
xct I fcnadd 
gs9bp, 

fadd sum 
store sum 



40.2 



mu Iby wf 9 



idx k 
sas (4 
jmp gs9l 

lio fcnadd+7 
Ipf 



load sum 
mulby hafwd 
jmp I gauss9 



54 , / 
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/phlnvr 

/inverse^of normal distribution, for statistical problems 

/Hasting s approximation 

/max error about 0.00004, which Is rather gross, but 

/I do not know of a better one at the moment. 

/uses flag 2 

/call by jdp phinvr with arg in AI, returns value In AI 
/note needs fin and fsr 

phlnvr, 

stf 2 
/check If < 0, out of range 

sma+sza 

jmp phnr7 

lac (400000 

c I i 

jmp i phinvr 
/check AI < 0.5 if so, use symmetry 
phnrf, dac phnrc 

sub (200200) /floating monotone in AI 

spa 

jmp phnrok 

elf 2 

lac phnrc 

negate 

fadd one 

dac phnrc 

/test if > 7, out of range 
sma^sza 
jmp phnrok 
lac (377777 
c II Icml 
jmp i phinvr 

phnrok, lac phnrc 

jdp fin 
mulby m2 
jdp fsr 
store phnrc 

/now the approx 

mulby phnrc+74 
fadd phnrc+72 
mulby phnrc 
fadd phnrc+70 
mulby phnrc 
fadd one 
store phnrc+76 
load phnrc 
mulby phnrc+6 
fadd phnrc + 4 
mulby phnrc 
fadd phnrc+2 
divby phnrc+76 
negate 
fadd phnrc 

szf 2 



Floating point package 



54.2 



negate 

jmp i phinvr 



phnrc, 
201240 
200315 
175251 
200667 

177301 
173653 




774353 
417430 
155426 
313144 
637356 
342312 
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/sort 

/ascending Shell sort 

/rorywrote and tested this on 7 feb 69 

/use jdp sortac with (A) = base address of vector to be sorted, 

/and (I) = number of them 

/saves X-register and flags/runs in lam, of course 



sortac, 





dio 


;nsort 


dio 


;msort 


sub 


(2 


dac 


;sortad 


TXA 




dac 


;savx 


20sort, 


lac msort 


sar 


7s 


dac 


msort 


sza 




jmp 


30sort 


Ixr 


savx 


jmp 


i sortac 



30sort, lac nsort 

sub msort 

dac ;ksort 

law 7 

dac ;j sort 
47sort, lac jsort 

dac ;isort 
49sort, lac isort 

add msort 

dac ; Isort 

/compare to see if need to swap 
lac sortad 
add Isort 
add Isort 
dac ;sort Id 
lac sortad 
add isort 
add Isort 
dac ; sort Id 
TAX 

/it would be faster to compare the upper and lower parts separately, 
/but more painful 

load i 

store so r tern 

Ixr sortld 

load i 

negate 

fadd sortem 

spq 

jmp 6 sort 

/swap 

50sort, Ixr sortid 

load i 
dimension sortem(2) 

store sortem 



Floating point package 



55.2 



Ixr 


sort Id 


loac 


i i 


Ixr 


sort Id 


store i 


Ixr 


sort Id 


loac 


1 sort em 


store I 


lac 


I sort 


sub 


msort 


dac 


I sort 


smaJ 


sza 


jmp 


49sort 


6 sort, 


idx 


sub 


ksort 


spq 




jmp 


47sort 


jmp 


2 sort 



jsort 



Floating point package 6 * 7 

/floating line display 

/property of rory thompson - written by charles landau 

/display a line from x7,y7 to x2,y2 

/limits of scope are _+7.0 

will be truncated to 7,0 



/all 


coord ina tes 


/call by 




/jdp 


fds 




/add 


ress 


of x7 


/address 


of y7 


/ditto xi 


I 


/ditto y2 


/control 


returns 


fds. 



TXA 






dac 


;savx 




law 


fdq 




dac 


fdp 




jmp 


fdp i 


fdq, 


dac 


dxT 




jdp 


fdp 




dac 


dyJ 




jdp 


fdp 




dac 


dx2 




jdp 


fdp 




dac 


dy2 




jdp 


dsp 




Ixr 


savx 




jmp 


i fds 


fdp, 









Ixr 


fds 




Ixr 


i 




idx 


fds 


fdr, 


lac 

spa 


i 


fix, 


cma 






sar 


Ss 




sub 


(400 




szm 






jmp 


f 7z 




add 


(43 




spa 






c la 






rar 


3s 




add 


(add f7y 




dap 


f7y 7 




and 


f7x 




ra I 


3s 




add 


(fa4 2 




dap 


fly 




lac 


i 




SXX 






lio 


i 




scl 


9s 


fly, 


xct 


• 




jmp 


• 




repeat 4, sc 




jmp 


i fdp 


fU, 


lac 


(377777 



here 



1377 makes scope limits 7/2, 401 ■* I, 



Ss 



Floating point package 56#2 



lio I 
sp! 

cma 

jmp I fjp 



Floating point package 



57.7 



/line dlsp lay 

^ d ? d , S ? Wlth coor J | nates of endpolnts in dx/ , dy7 , dx2, dy2 
/dx7,dy7,dx2,dy2 are saved 
/xr is not saved 



dx7, 

dy7, 

dx2, a 

dy2, 



dx3, 
dsp, 



csp / 



Jy3, 



dn / 



el n, 



/x coordinate of endpolnt 
/y coordinate 



dx2 
dx7 
7 
7 



/cou Id be modif ed 
/x2 and x7 in AI 







lac 

lio 

sar 

si r 

AMIX 

lac dy2 

lio dy7 

sar 7 

si r 7 

AMIAI 

dac dy3 

spa 

cma 

TXIX> 

AMI I | /| x2-x7| +1 y 2-y7 

A+II 

c la 



to enter with 



scr 

scr 

dio 

TXA 

CIX 

jmp 

mu I 

div 



sa I 

dac 

lac 

mu I 

di v 



sa I 

dac 

lac 

lio 

iot 



2s 

6s 
dn 



/increase for wider spacing 



cin-2 
(7 

dn 

7 

dx3 

dy3 

(7 

dn 

7 

dy3 

dx7 

dy7 

207 



SXX<M 

jmp i 

swp 

add 

swp 

add 

jmp 



dsp 

dy3 

dx3 
ci n 
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t P .K Sr ?n D T?AM n »?! rate s y stens of differential equations is available 
in the FORTRAN library. It uses a Kutta-Merson scheme with automatic 
step size adjustment, 

ca I I ktmbeg(array, deriv, n, err, I i m. In it, max) 
sets up the arrays and integration parameters. 

Al I are rea I except n. 

Array is a real array containing the variables. ArrayO ) is the 
independent variable. 

Deriv Is a real array containing the derivatives of the above array. 
DerivO) must contain 7.0. The calling program is responsible 
for this. 

N is the number of variables in array and deriv. 

Err Isthe maximum allowable error In any step of the integration. 
If it is exceeded, the program will do the step over with a smaller 
step size. 

Lim is the limit of Integration for the independent variable. The program 
will stop when the variable has increased by this amount, regardless 
of its initial value. 

Init is the initial step size. It may be adjusted later if necessary. 

Max is the maximum step size. When adjusting the step size, It will 
never be permitted to exceed max. 

ca I I ktm 

This begins the integration. The integrator will periodically (about 
four times per step) ask for the derivatives of the variables in 
array. The first time it will do so by returning from ktm. The 
calling program must then compute the derivatives of the current 
values of array and store same in deriv. It then returns to the 
integrator by 

call ktmret(70s) 

where 70s is a statement number followed by the letter s (an 

abnormal subroutine return N argument). If the Integrator wants 
more values of the derivatives, it will return via this argument. 
The calling program must compute the derivatives again and call 
ktmret again. When ktmret makes a normal exit, the Integration 
is complete, and the results are in array. It Is then possible 
to call ktm again, and integrate further by a distance of lim 
In the independent variable. 



Floating point package 



60.1 



examp le 



integration at zero 

right constant of integration 



dimension y(2),d(2) 

call ktmbeg(y,d,2,7e-7,. J,.0 7,7.) 
c initial step size = .07, stop every 

d(7)=7.0 

y(7) = 0.0 
c to begin the 

y(2)=0.0 
c to give it the 
10 call ktm 
10 d(2)=7.0/(7,0+y(7)**2) 

ca 1 1 ktmret(20s) 
c dy/dx = 7/(7+xt2) 

cal I fop(y(7)) 

call tab 

cal I fop(y(2)) 

ca 1 1 ret 
c this makes a table 

go to 10 



of arctangents in steps of 0.1 



Floating point package 

/Kutta-Merson 

/before calling, define Initial values y(0) to y(n) 

/where y(0) is the independent variable, 

/define fin, the end of the domain, 

/define n, the number of equations, 

/define err, the approx max allowable error, 

/not to mention h, the initial step size 

/(which will be adjusted) 

/or him, the maximum step size 

/fen is a function that returns the derivatives 

/of y(n) in kfn(n). i.e., kfn(0)=7 

/int, the address of a location to handle the intermediate 

/results ( none, If zero). then jdp ktm 

siz=20. 

nm=2*si z 

dimension y0(nm), f 0( nm), f / (nm) 

dimension err(2 ), hc(2 ), h(2 ) 

dimension hc6 (2 ), hc2( 2 ) ,epv(2 ) 

dimension hcS(2 ), f in( 2 ) 

dimension hc3(2 ),h lm(2 ) 

n, /number of equations 
int, /address 

ktm, /use jdp 

iam /now can run In Index mode 



Floating point package 6 Z * 7 

/p2 of ktm 

I io n 
Law si z 

AMI> = 

bpt /redefine siz and assemble again 
si I 7s 
dio ;2n2 

load h 
store he 

km0, lac i y 
dac i y0 
SXAX 
sas 2n2 

jmp km0 

load err 
divby ten 
store epv 

km/, load he /calculate these factors but once 
divby thr 
s to re he 3 
sub (400 
store hc6 
load he 
sub (4 00 
store hc2 
sub (1000 
store hcS 

/test if done, if fin - y0 < err 
kmc, load y0 

cma+cmi *c If 6 

fadd fin 

store tern 

negate 

fadd err 

sma 

jmp i ktm /exit 

/if too close to end ( fin - y0 < he), change he 
load tern 

fadd epv /to watch out for equal giving minus zero 

negate 

fadd he 

spa 

jmp kms /he okay 

load tern 

store he /change he 

jmp km7 /get other constants 



Floating point package 53,7 

/p 3 of ktm 
/proceed with step 
kms, jdp fen 

2X 
km2, load i kfn / Lkf n] $ [y'(t)] ■* tfQl 

store ? fO /[kfn]h/3 + [y0]-*[y] 

mulby hc3 
store tern 
load i y0 
fadd tern 
store ? y 
SXX 
SXAX 
sas 2n2 
j mp km2 



Floating point package ' 

/ p 4 of ktm 

jdp fen 

£X 
km3, load i kfn /([kfn] + [f0] )h/6 + [yfl] ■* [y] 

store tern /[y] $ [y(t+h/3)] 

load i f0 
fadd tern 
mulby hc6 
store tern 
load I y0 
fadd tern 
store I y 
SXX 
SXAX 
sas 2n2 
jmp km 3 
jdp fen 

2X 
km4, load i kfn /[kfn] $ Cy'(t+h/3)] - [f7] 

store i f7 /([kfn33 + [f0])h/S + [y0] - [ y] 

mulby thr /[y] $ [y(t+h/2)3 

store tem 

load i f0 

fadd tem 

mulby hcS 

store tem 

load i y0 

fadd tem 

store i y 

SXX 

SXAX 

sas Zn2 

jmp km4 

jdp fen 

2X 
km5, load f kfn /[kfn] $ [y'(t+h/2)] -> [f2] 

/([kfn] 4 - [f/]3 + [f0])h/2 + [y0] -* [y] 

mulby for / [y] $ [y(t+h)] 

store tem 

load I f0 

fadd tem 

store i f0 /save 4kfn + f0 In f0 for next step 

store tem 
load I n 
mulby thr 
negate 
fadd tem 
mulby hc2 
store tem 
load I y0 
fadd tem 
store I y 
SXX 
SXAX 
sas 2n2 
jmp km5 



Floating point package 

/ page 5 of ktm 

jdp fen 

2X 
km6, load i kfn 
store tern 
load i fO 
fadd tern 
mulby hc6 
store tern 
load I y0 
fadd tern 
store I f7 
store tem 
load i y 
negate 
fadd tem 
spa 

negate 
sub (1000 
store tem 
negate 
fadd err 
sma 
jmp km7 



6 5.7 



/([kfn] + [f0] + lf214)h/6 + [y0] 
/if 11 $ [y(t+h)3 



[f/] 



/test if really want this step 
/test if should change step size 

/dlvby 4 

/given error bound 



6 6 1 
Floating point package 



/ page 6 of ktm 

law i 400 /error too big 

adm he 

2X 

km9, lac i y /reset for first call of fen again 

dac i y 
SXAX 
sas Zn2 

jmp km 9 

j m p km 7 

/ here if error acceptable 
km7, load tern 
negate 

fadd epv /check if small enough so can increase steo size 

spa 

stf 6 /dont increase size 

SXX 

SXAX 
sas 2n2 
jmp km6 

=X /accept the step 

km£, lac i f7 

dac i y 

dac i y0 
SXAX 

sas 2n2 

jmp km S 

lac int /vent intermediate to outside 

szm 

jdp i int 

kmz, szf 6 

jmp kmc 

load he 

add (4 00 /*Z 

negate 

fadd him /for guard against too large a time step 

spa 

jmp kmc 

Law 4 00 

adm he 

jmp km) 

/end of ktm 

/remember to define int If use it. Alternatively, can set fin at 

/various desired steps and output then. 



10 tape 
bs=l0 

100/ hit 
hit 

law 500 
lio (30001 
ivk 16 
hit 

rdd, lio (dir 
law 100 
jdp rdb 
jmp cr-1 

consta 



cr-1, 
cr, 



tbl, 



1st, 
tb2, 



err 



iam 
law 15 
jdp atof 
ivk 300 
jdp rd6 
TIXP 
jmp err 
dac 1st 
sad 1 tbl 
xct 1 tb2 
SXX 
jmp .-3 



420000 
440000 
450000 

460000 
560000 





jmp 
jmp 
jmp 
jmp 
jmp 
jmp 
jmp 



or 

Id 

del 

edl 

fll 

err 

err 



law 74 
law 21 
law 72 

jmp cr 



/null command 

/b - return to id 

/d - delete file 

/e - edit file 

/f - file file 

/p - print directory 

/error 



/print not implemented 



ivk 100 
ivk 100 
ivk 100 



rdb, 



wtb, 



trb, 








jdp 


trb 


ivk 


16 


hit 




jmp 


i rdb 







jdp 


trb 


ior 


(add 


ivk 


16 


hit 




jmp 


1 wtb 







rar 


Is 


sma 




cma 




and 


(777 


swp 




jmp 


i trb 



/read sixbit name into AC-IO and nml-nm2 
/call with jdp rd6 
rd6, 

elf 7 

dzm nml 

dzm nra2 

ZX 
rdl, ivk 400 

jdp ftoa 

sub (40 

spq 

jmp dun 

and (77 

cli 

aam 

jmp .+1 
repeat 5jrcl 6s 

A-»IA| 

ZAI| 

SXX 

adm nml /yes, I know 

lai 

adm nm2 

jmp rdl 
dun, sad (10-40 

jmp cr 

lac nml 

lio nm2 

jmp i rd6 



gfn, 









szf 


1 




jmp 


err 




jdp 


rd6 




sza 






szf 


1 




jmp 


err 




dac 


nam! 




dio 


nam2 




Jdp 


rd6 




szaVszf i 1 




jmp 


err 




law 






dap 


fwd 




jmp 


1 gfn 


ffn, 









dzm 


Fnl 




ZX 




ffl, 


Idx 


fnl 




lac 


naml 




sas 


1 dlr 




jmp 


gfi 




lac 


nam2 




sas 


1 dir+1 




jmp 


gfi 




lac 


nml 




sas 


1 dlr+2 




jmp 


gfi 




lac 


nm2 




sad 


1 dlr+3 




jmp 


fff 


gfi, 


law 


4 




X+AXA 




sas 


(4x27 




jmp 


ffl 




jmp 


1 ffn 


fff, 


idx 


ffn 




jmp 


i ffn 



/return to id 

id, 

wd, skp i 600 
jmp id2 
law 100 
lio (buf 
jdp rdb 
lac dir+377 
sas buf +377 
-i 

law 100 
lio (dir 
jdp wtb 

id2, law i 
ivk 16 
nop 
law 16 
rata 204 
dsm 



fil, 



rce, 
rce+3> 



nmf 9 



jdp 


gfn 


law 


2 


jdp 


rset 


jdp 


ffn 


jmp 


.+2 


jmp 


err 


dzm 


fnl 


szs 


20 


stf 


2 


claVclf 1 


dap 


wd 


jdp 


rev 


ZX 




dzm 


Fn2 


idx 


fn2 


lac 


1 dir 


ior 


i dir+1 


ior 


1 dir+2 


lor 


i dlr+3 


sza 


1 


jmp 


nmf 


law 


4 


A+XXA 


sas 


(4x27 


jmp 


rce+3 


jmp 


err 


lac 


naml 


dac 


1 dir 


lac 


nam2 


dac 


1 dlr+1 


lac 


nml 


dac 


1 dir+2 


lac 


nm2 


dac 


i dlr+3 


dzm 


naml 


dzm 


nam2 


dzm 


nml 


lac 


fn2 


dac 


nm2 


ral 


Is 


TAX 




law 


i 1 


and 


i dir+2x56-l 


dac 


i dir+2x56-l 


law 


1 1 


and 


i dir+3x56-l 


dac 


i dir+3x56-l 


law 


1 


lor 


i dir+4x56-l 


dac 


i dir+4x56-l 


jdp 


fnb 


jmp 


err 


jmp 


flll-1 



fnfr, 


law 1 bs-1 




xct fwd 




adm blk 




jdp fnb 




jmp rce 


fill-1, 


dzm 7"pt 


fill, 


law 12 




szs 1 30 




szf 1 1 




jmp .+2 




jmp f!13 




jdp readch 




jdp ftoa 


fil2, 


jdp put 




idx 7pt 




sas (200x5 




jmp fill 




lac blk 




llo (buf 




jdp wtb 




law 3 




sas to 




jmp fnf r 




jmp cr 


fil3, 


elf 1 




jmp fll2 



edi, 


jdp 


gfn 




law 


2 




jdp 


wset 


te, 


jdp 


ffn 




jmp 


err 




jdp 


rev 




dzm 


naml 




dzm 


nam2 




dzm 


nml 




lac 


fnl 




dac 


fn2 




dac 


nm2 


edb, 


jdp 


fnb 




jmp 


te 




lac 


blk 




llo 


(buf 




jdp 


rdb 




dzm 


7pt 


edil, 


jdp 


get 




jdp 


atof 




jdp 


wrltec 




idx 


7pt 




sas 


(200x5 




jmp 


edll 




jmp 


edb 



del, 



dee, 



rev, 



jdp 


gfn 


law 





dap 


wd 


dzm 


fn2 


jdp 


ffn 


jmp 


cr 


dzm 


1 dlr 


dzm 


i dlr+1 


dzm 


1 dir+2 


dzm 


1 dlr+3 


dzm 


naml 


dzm 


nam2 


dzm 


nml 


lac 


fnl 


dac 


nm2 


jdp 


rev 


jdp 


fnb 


jmp 


dee 


jmp 


.-2 







law 


cma-opr 


xor 


fwd 


dac 


fwd 


law 


400 


xct 


fwd 


add 


(400 


dac 


blk 


jmp 


1 rev 



fnb, 
fwd, 



5st, 



5ut, 








law 


i 1 


opr 




a dm 


blk 


lio 


(1000 


szm 




AMI< 


jmp 


1 fnb 


CSC 




mul 


i? 


div 


hit 




dac 


to 


TAX 




lac 


1 dlr+4x27 


lxr 


i dir+4x27+l 


X->IX 


xct 


i 5st 


dac 


tl 


and 


(37 


sas 


fnl 


jmp 


fnb+1 


law 


i 37 


and 


tl 


lor 


fn2 


xct 


1 5ut 


lxr 


to 


dac 


1 dir+4x27 


dlo 


1 dlr+4x27+l 


Idx 


fnb 


jmp 


i fnb 


ral 


5s 


ral 


7s 


rar 


8s 


rar 


6s 


rar 


3s 


rar 


Is 


rcl 


2s 


rar 


5s 


rar 


7s 


ral 


8s 


ral 


6s 


ral 


3s 


ral 


Is 


rcr 


2s 



fcs, sub (72 
ral 5s 

dac cas 
law i 2 
adm .+2 
jmp i .+1 

ftoa, 

laiVclf 1 
sas (72 
sad (74 
jmp fcs 
ior cas 
add (tab 
dap .+1 
lc 

and (177 
sad (15 
stf 1 
szf 2 

jmp i ftoa 
sub (140 
spa 

add (40 
add (lOO 
jmp i ftoa 

atof, 

sad (12 
jmp If 
sad (15 
jmp .cr 
elf 6 

af, TAAX 

sub (100 
szs 50 
TAA> 
jmp afl 
sub (32 
TA> 
law 40 
A+XX 

afl, lac 1 tab 
ell 

rcl 9s 
law 100 
AAIA 
sad cas 
jmp 1 atof 
dac cas 
rar 5 s 
add (72 
swp 
aam 

xct atof 
11a 
jmp 1 atof 

If, szf 1 6 
jmp af 



cr 



idx at of 

elf 6 

jmp i at of 

szs i 40 
stf 6 
jmp af 



8 
get, 



put, 



Pg2, 



Pg3, 
sht, 



ust, 








elf 


4 


jmp 


Pgl 







dac 


to 


stf 


4 


lac 


7pt 


mul 


(146315 


dac 

AMI] 


tl 

[ 


TAX 
rcl 


3s 


sza 


1 


law 


7 


llo 


1 buf+1 


lxr 


i buf 


X^AX 


xct 


1 sht 


szf 


1 4 


jmp 


Pg3 


and 


(-177 


lor 


to 


xct 


1 ust 


lxr 


tl 


dac 


1 buf 


dio 


1 buf +1 


jmp 


i put 


and 


(177 


jmp 


1 get 


hit 




rcl 


3s 


hit 




rar 


Is 


rar 


4s 


hit 




rar 


8s 


ral 


7s 


hit 




rcr 


3s 


hit 




ral 


Is 


ral 


4s 


hit 




ral 


8s 


rar 


7s 



define z f,a 




f xlOOO a 


termin 




tab, z 


76, 4o 


z 


76,61 


z 


76,62 


z 


14,63 


z 


76,64 


z 


76,65 


z 


76,66 


z 


76,67 


z 


75,70 


z 


36,71 


z 


77,0 


z 


76,14 


z 


13,3 


z 


77,0 


z 


34,0 


z 


35,0 


z 


76,60 


z 


76,57 


z 


76,163 


z 


76,164 


z 


76,165 


z 


76,166 


z 


76,167 


z 


76,170 


z 


76,171 


z 


76,172 


z 


76,0 


z 


76,54 


z 


76,17 


z 


76,16 


z 


76,11 


z 


76,0 


z 


0,72 


z 


105,152 


z 


101,153 


■ ^ z 


103,154 


z 


103,155 


.- *>z 


104,156 


z 


106,157 


z 


102,160 


z 


57,161 


z 


55,162 


z 


173,0 


z 


154,0 


z 


33,55 


z 


54,51 


z 


73,73 


z 


21,50 


z 


20,0 


z 


1,141 


z 


2,142 


z 


3,143 


z 


4,144 


z 


5,145 


z 


6,146 


z 


7,147 



~ - ~> 



10,150 

11,151 

4o,o 

56,56 

107,0 

133,10 

110,0 

121,15 

120,40 
161,42 

162,47 
163,44 
164,45 
165,41 
166,46 

167,74 
170,76 
171,136 
l4l,o 

142, Oil 

143,03 
144,0 
145,0 
z 146,0 
z 147,100 

150,77 

151,123 

122,124 

123,125 

124,126 

125,127 

z 126,130 
z 127,131 
130,132 
131,0 
157,75 
156,17 
155,16 
111,11 
l4o,o 

102,137 
61,112 

62,113 

63,114 

64,115 

65,116 

66,117 

67,120 

70,121 

71,122 

4l,o 

42,0 

43,53 

44,135 

45,134< 

46,133 
47,0 
50,101 
z 51,102 
z 22,103 



/100 



z 23,104 

z 24,105 
z 25,106 
z 26,107 
z 27,110 
z 30,111 
z 31,0 
z 157,52 
z 156,0 
z 155,10 
z 103,0 
z 76,15 



rbl=400 

dimension rbf(rbl) 

rset, 

dap rsr 

dac ac 

dio io 

rar 6 s 

add (400-rbl 

dac ddp__ 

dzm rpc 

lac (isp rpc 

dac rfg /flag for end of file 

law rbf+rbl-1 /last word of buffer 

dap inp 

law rbf 

mta 

lio (3^0 

law 40 
rsr, ivk . 

hit 

lac ac 

rar 6 s 

add rbf 32 

sub (20000 

dac dne 

lio io 

lac ac 

jmp i rset 



2 








readch, 











dac 


ac 




rfg, 


isp 


rpc /becomes 


jmp rp2 after 




jmp 


inp 






law 


i 3 






dac 


rpc 






idx 


inp 






dap 


inq 






sas 


(lio rbf+rbl 






jmp 


9k 






lw rbl 






a dm 


ddp 






11a 








ral 


6s 






dap 


rdr 






law 


rbf 






dap 


inp 






dap 


inq 






mta 








law 


rbl 




rdr, 


Ivk 
hit 

ZAP 






9k, 


sub 


(lio rbf 






add 


ddp 






sad 


dne 






jmp 


rp3 




inp. 


lio 
cla 




: 




rcl 


6s 




inq,. 


dlo 








sas 


(77 

(13 






sad 






dzm 


rpc 




rpl, 


lia 








lac 


ac 






jmp 


i readch 




rp3, 


lac 


(jmp rp2 






dac 


rfg 




rp2, 


lw 14 






jmp 


rpl 





end of file 



ofi=4oo 

dimension ofb(ofl) 

wset, 

dac ac 
dap wrdr 
mta 300 
nop 
sni 
bpt 

rar 6 s __ 
dac wfld 
add _(400-of 1 
dac odp 
law i 3 
dac cp 
law ofto 
dap ofp 
dap ofq 
lac ac 
jmp i wset 



writec, 



wo5> 



wrd,- 

ofp-1, 
ofp, 



ofq, 



wox 



ou3, 



dio 


95p 


dac 


ac 


law 


77 


AAIA 


sad 


(14 


jmp 


clean 


dac 


?6p 


lac 


(lio ofb+ofl 


sas 


ofp 


jmp 


ofp-1 


law 


ofl 


a dm 


odp 


11a 




and 


(770000 


sas 


odp 


jmp 


wo5 


ral 


6s 


mta 


300 


nop 




snl 




bpt 




jmp 


.+2 


ral 


6s 


dap 


wrd 


law 


ofb 


dap 


ofp 


dap 


ofq 


mta 




law 


ofl 20 


ivk 




hit 




lac 


96p 


lio 




rar 


6s 


rcl 


6s 


dlo 




isp 


cp 


jmp 


ou3 


law 


i 3 


dac 


cp 


Idx 


ofp 


dap 


ofq 


lac 


ac 


lio 


95p 


jmp 


1 wrltec 


lac 


96p 


sas 


111 

(13 


sad 


jmp 


ofp 


jmp 


wox 



/write out buffer 



clean, law 13 
lia 



sas 96p 

jdp wrltec 

lac ofp 

add odp 

sub wfld 

add (20000+ofl-[llo ofb] 

dac b77 

lw ofb+ofl 

dap ofp 

jdp writec /write out last buffer 



law b77-32 

mta 

lio (3^0 

law 60 
wrdr, ivk . 

hit 

jmp cr 
b77, 

20400 

constants 

epeat ifm 2100-. ,printx /lose/77// 

variables 



5000/ 



ini, 


ZAP 


in, 


lw 600 




dap inn 




lac (add 100 
lio (40000 






ivk 16 




hit 


inn, 


law 




dap wd 




sza 




jmp rdd 




jmp cr-1 


constant 





define w a,t>,c,d,e,f ,g 
ax20000 bx400 cxlO [dA34J>4 
[dA3]x200000 ex4000 fxlOO gx2 
Vermin 

6000/ 

dir, repeat 134,0 

w 36,36,36,36,36,36,36 

repeat 2x10,0 

w 33,0,0,0,0,0,0 

repeat 2x77,0 

repeat 10, w 36,36,36,36,36,36,36 

-0 

675756 

repeat ifn ,-6400,printx /loses/77// 

buf=7000 

start ini 



